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The dependence of inter-particle correlations on the orientation of particle relative-momentum can 
yield unique information on the space-time features of emission in reactions with multiparticle final 
states. In the present paper, the benefits of a representation and analysis of the three-dimensional 
C ■ correlation information in terms of surface spherical harmonics is presented. The harmonics include 

the standard complex tesseral harmonics and the real cartesian harmonics. Mathematical properties 
of the lesser-known cartesian harmonics are illuminated. The physical content of different angular 
harmonic components in a correlation is described. The resolving power of different final-state effects 
with regarding to determining angular features of emission regions is investigated. The considered 
final-state effects include identity interference and strong and Coulomb interactions. The correlation 
analysis in terms of spherical harmonics is illustrated with the cases of gaussian and blast-wave 

■ sources for proton-charged meson and baryon-baryon pairs. 

<n : 

(N) ■ I. INTRODUCTION 

> : 

' Measurements of low relative- velocity correlations yield access to the size and shape of relative emission sources in 
reactions with multiparticle final states [ll, 0> S 01 ■ This information plays an especially important role in in heavy- 

■ ion collisions. The duration of emission, which is linked to the nuclear equation of state, can be also inferred from 
t-H \ those measurements. Perhaps the most startling finding, during the first years of the Relativistic Heavy Ion Collider 

(RHIC) program, was the disagreement between the features of experimentally determined emission sources [H, H, 
and predictions of hydrodynamic models [1, [ij [13, HH incorporating an anticipated phase transition, which would yield 
a long emission duration. A long emission duration would have resulted in an elongated shape (i? ut /-Rside >> 1 in 
the usual RHIC nomenclature) of the emission source for particles of a specific momentum. In the literature, the most 
exhaustively studied correlation has been that produced by interference for identical mesons. For pure interference, 
the correlation expressed as a function of the relative particle momentum is related to the relative emission source 
through a simple Fourier transformation, principally yielding direct access to the source shape. In the measured 
charged-meson correlations, the Coulomb effects competing with interference had been normally compensated for 
with cumbersome corrections. However, as we demonstrate, the correlations induced by the Coulomb and further by 
the strong interactions can also provide information on the source shapes. For any type of correlation, the analysis 
is facilitated by the use of spherical harmonics, either the standard complex tesseral harmonics or the real cartesian 
harmonics [l2T ] . The harmonics allow for a straightforward faithful representation of the three-dimensional correlation 
data, and also better facilitate isolating and focusing on specific physical properties of the source. Different source 
parameterizations will be employed for illustration of our discussions. For example, the I — 2 harmonic expansion 
coefficients, characterizing quadrupole deformation, provide information on the ratio of radii i?out/-Rsidc in the common 
gaussian source representation. 

The correlation, associated with the final-state effects at low relative velocities within a subsystem of the particles 
a and b, can be linked to a source function S through the relation: 
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■K(P,q) = J d z r\<j) { -\q,r)\ 2 S{P.r). (1) 



Here, 1Z(P, q) is the deviation of the l.h.s. correlation function from unity, P is the total momentum of the pair 
and q is the relative momentum measured in the two-particle rest frame. The factor |0( _ -'| 2 represents the relative 
wavefunction squared, in the two-particle rest frame, with asymptotic condition imposed on the outgoing wave. The 
square of the wavefunction is averaged over spins. The function S{P,r) represents the probability density that the 
particles a and b, moving with the same velocity V corresponding to the total momentum P, are emitted from the 
reacting system a distance r apart in their rest frame. Since at large q, the emission in multiparticle final states 
becomes uncorrelated, i.e. 1Z — > while l^ - -*! 2 averages to 1 as a function of r, it follows that iS is normalized to 1, 
/ d 3 r S(P,r) = 1. A detailed derivation of Eq. ([T]) can be found in [l3j] (see also [bj]). 
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The probability density of emission at the relative distance may be expressed in terms of single-particle probability 
densities s: 

S(P r) - / d 4 z Q d 4 x b s a (m a P/(m a + m b ),x a ) s b (m b P/(m a + m b ),x b ) 5{r -r a -r b ) 
J d A x a d 4 x b s a (m a P/(m a + m b ), x a ) s b (m b P/(m a + m b ), x b ) 

where s a (p, x) is the density for emitting particles of type a with momentum p from space-time point x. The spatial 
arguments are for the frame moving at the center-of-mass velocity V . According to Eq. |T]), the correlation data access 
the probability distribution in the relative rest-frame separation only. Information on emission duration can only be 
inferred from the influence of that duration on the relative spatial distribution [IH [Til [ItJ • Similarly, single-particle 
distributions s can be deduced only indirectly from data |T]) because of the folding in ((21) with an integration over the 
average particle position. For a set of investigated particles, the constraining of single-particle distributions may be 
enhanced by an analysis of correlations ([T]) of all possible particle pairs. 

In the following, we concentrate on what might be directly determined from data, i.e., the determination of S(P,r) 
given measurement of 7?.(P, q) . Since each individual value of P nominally represents a different problem, we suppress 
the P arguments in the functions. As access to S is conditioned on any deviations of \(f)\ 2 from unity, we subtract 
unity from both sides of |T]) arriving at 

K(q) = j d 3 r [\4> { -\q,r)\ 2 - l] S{r) = J &\ K{q,r) S{r) , (3) 

where we have utilized the normalization of S(P,r). Here, both q and r are evaluated in a frame where P = 0. It 
is seen that 1Z is related to S through an integral transform where \4> 2 (q,r)\ 2 — 1 plays the role of the kernel JC of 
transformation. The determination of S{r) from lZ(q) amounts then to the inversion of the integral transform, often 
referred to as imaging [H, [H H3, HI [H, Q . 

Imaging has already been extensively performed for angle-averaged correlation functions, lZ(q), which can provide 
angle-averaged distributions S(r). It is the goal of this paper to describe how one might extract the angular infor- 
mation in <S(r) by decomposing lZ(q) in terms of surface-spherical tesseral and cartesian harmonics. A parallel effort, 
employing tesseral harmonics in the case of pure interference, is described in 24]. 

The wavefunction squared \cf>(q,r)\ 2 depends, after summation over spins, only on the magnitudes of q and r and 
on the angle 9 qr between q and r. Thus, the square is invariant under rotations. As a consequence, the coefficients 
of expansion in terms of spherical harmonics, for the correlation and source functions, are directly related to each 
other [13, 



n lm {q) = Air J dr r 2 K t {q, r) S £m (r) . (4) 
Here, the expansion coefficients and expansion for the correlation function are defined by the relations: 
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Tl tm {q) = -7!= / dn q Y hn (Q)n(q), (5) 



4ir 

with analogous expressions for the source function. The partial kernel ICi encodes the information from the wavefunc- 
tion, 

fce{q,r) = ^ J dcos8 qr (q,r, cos 9 qr )\ 2 - 1 P e (cos9 qr ) . (6) 

According to the above, an (£, m) coefficient of the correlation function is tied to the coefficient of the source function 
with the same (£,m). Thus, the source restoration can be reduced to a series of one-dimensional inversions of the 
integral transformations for individual (£,m), employing kernels for specific £. This sequence of one-dimensional 
inversions is computationally simpler than a single three-dimensional inversion. The angular resolution of data would 
determine how far the restoration can proceed in £. 

One problem with the above is the emergence of complex coefficients for the real-valued correlation and source 
functions, lacking direct geometric interpretation. The last issue is circumvented by introducin g an alternate real- 
valued basis in the space of spherical angles, surface-spherical cartesian harmonics A^i) [l2l. l25l. I2H] . 
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Cartesian harmonics Ag are linear combinations of Yi m with different m but the same £ = £ x + £ y + £ z . In terms of 
the directional unit vector n pointing in the O direction, the cartesian harmonics are explicitly given by [l2| 



l\ m (2£-2m-l)!! £ x ! 



<V V 2 J (2£ - 1)!! (4 -2m x )\m x \ 

m 

0<m,i<t.i/2 

/I /I 

(£ y -2my)\m y \ (£ z - 2m z )\ m z l v [) 

where m — m x + m y + m z and (—1)!! = 1. The leading term of „4^*is n l * n y n z ", while the subsequent terms ensure 
that r l Ag satisfies the Laplace equation and, thus, is a combination of the Y( m of different m. 
The expansion in terms of cartesian harmonics are defined with the relations 

= £ 777777 KtMO*) = £ 777777 n r& & & ' ( 8 ) 
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In many situations the explicit form of .A^* is not needed. Just the existence of the cartesian harmonics and their 
properties justify representing a function of the spherical angle in the second series in ([8]). 

Since Ag is a combination of Yi m with £ — £ x + £ y + £ z , the coefficients of expansion in terms of cartesian harmonics 
for the source and correlation function are related to each other in the same way as the coefficients of expansion in 
terms of tesseral harmonics in Eq. ([4]), 

K £ iq) =4irjdr r 2 JC e (q, r) S^r) . (10) 

In the next section, the structure of the kernel JCg in Eqs. ^ and (fit))) is investigated for the cases of pure identical- 
particle interference and for strong and Coulomb interactions within a pair. We find that all three classes of final-state 
effects provide significant resolving power for different £. Section IlIII provides a detailed discussion of the properties 
of cartesian harmonics. The subsequent section presents a discussion of the utility of higher £- moments and their 
relation to specific geometric features of the source. In particular, we show that the £ = 1 moments can reveal an 
offset between the probability clouds for two different species, that the £ = 2 moments can reveal the magnitude and 
orientation of axes in the ellipsoidal approximation for a source, while the £ — 3 moments can reveal the "boomerang" 
nature of a source. In Sec. |Vl the cases of gaussian and blast wave sources are studied, as examples, in quantitative 
detail. In the final section, we summarize our results and discuss the prospects and challenges expected in analyses 
employing the harmonics. 



II. CORRELATION KERNELS K t 



The ability to extract source S anisotropics depends on the kernels JCi{q,r) defined in Eq. Data on 1Z are 
obtained at a certain resolution in q and suffer from errors. The theoretical relation |T]) involves approximations [l3l[2^ |. 
in assuming that the relative momentum q is small compared to the scales characterizing the rest of the system and 
in ignoring final-state effects other than those between the two studied particles moving slow relative to each other. If 
the kernels ICi drop rapidly as £ increases, for typical q- values and for values of r characteristic to the source function, 
there will be no chance to determine details of the shape of S(r). In the following we analyze the situation for different 
types of final-state interactions at different £, q and r. 

In the process of source restoration, the values of IZj? would be known at discrete values of q. In parallel, the source 
Sjr(r) might be discretized in r or decomposed in some basis. The source restoration would then amount to the inversion 
of a matrix out of K,i{q,r), cf. pll. [28|. An inability to extract source features might be signalled by the proximity 
of matrix eigenvalues to zero, resulting in an instability of the inversion. Such a situation could be encountered 
when extracting details of in an r-region where ICi lacks resolving power. Practical experience in matrix-inversion 
strategy, in analyzing data, has been gained for £ = [18|, [l9|, |2CJ, [2l|, |22|, 23] . Importantly, this strategy has revealed 
non-Gaussian features in Se=o(r), which were especially significant for proton sources in intermediate energy collisions. 
More common in the literature than the source discretization, has been a parameterization of the source which is 
then fit to 1Z. When parameterizing the source in three-dimensions, the benefit regarding angular moments can be 
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in understanding a systematic of the moments, TZt and S/, as a function of £. With regard to the fit strategy, the 
features of ICe (g, r) determine what parameters of S may be potentially constrained by data. 

In the first subsection we discuss kernels for the case of pure identity interference between spin-0 bosons. The two 
subsequent subsections are dedicated to the Coulomb and the strong interactions. Before launching into a detailed 
discussion of specific classes of interactions, we emphasize a general observation that a three-dimensional source Sg(r), 
expanded in a Taylor series around r = 0, behaves as r for small r. Consequently, there is not much to be discerned 
about the source at small r for £ > 1, particularly for large £, and the structure of Kg>x at small r will not be 
overly important. On the other hand, much could be learned about the angular structure of the source at large r. 
However, when the large-r region gets mapped onto small g, the ability to gain angular information becomes limited 
by the experimental resolution capabilities. Thus, in practice, the most important features of K,e>i(q,r) are those 
at moderate q combined with moderate to high r. The naive expectation regarding strong interactions may be that 
when an s-wave dominates the interaction, there is no possibility to discern source shapes. We shall see that this is 
not the case. 



A. Kernels for Spin-0 Boson Identity-Interference 

For non-interacting identical spin-0 bosons, the squared relative wave function is 

\(f>(q, r, cos 9 qr )\ 2 = 1 + cos (2qr cos9 qr ) , (11) 

where the relative momentum convention is that q stands for the momentum of one of the particles in the center-of- 
mass frame. In this case, the three-dimensional transformation ([3]) becomes the cosine Fourier transform [23$ . Since 
the cosine term in (TTTj) may be represented as the real part of e 2iq r and the latter may be expanded in Legendre 
polynomials and spherical Bessel functions, the /C^-transformations (j!} and (JTUJ1 . with (J6]), emerge as spherical-Bessel 
transforms (see also [Hj]) : 



^-{(-^^•ZZtt: c* 

The properties of ji(2qr), as a function of qr, are as follows. At low qr, the function starts out in the power form, 
je(2qr) w (2qr) e /£\\, which is followed, for I > 2, by a maximum at qr ~ | + \ and by a subsequent alternation at the 
approximate period of tt. This has the following consequences for the ability to discern the source angular components 
Si. If R is a scale characterizing Si, Se will have little impact on IZi if, at £ > 2, qR falls into the hole in je at low 
argument values. Likewise, Se will have little impact on IZe , if many oscillations of je are averaged over in the spatial 
integration over the source. The maximal impact is expected if qR falls into the region of the first maximum in je 
and, thus, most information on Sg may be expected in He at q ~ {£ + ^)/(2R). Otherwise, the coefficients IZe are 
expected to behave as q for low q and vanish at large q. 

As characteristic for Fourier transforms, q and r enter the kernels as a product, and the low-r region of the source 
gets predominantly mapped onto the high-q region of the correlation function. Conversely, the high-r region of the 
source function mainly affects the \ow-q region of the correlation function. The inverse to the cosine Fourier transform 
is again the cosine Fourier transform and the inverse to the spherical Bessel transform is the spherical Bessel transform. 
Given the Bessel function completeness relation, 

dqq 2 je(2qr)je(2qr') = I-S(r-r'), (13) 
16 



we can introduce an inverse kernel 



and obtain 



IC7\r,q) = { ^U^qr) , for even £ , (w) 
L , for odd £ , 



Sfr) = 4tt / &qtfK7\r,q)TlM . (15) 
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FIG. 1: Kernel, IC(q,r,cos8 q 
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1, for classical repulsive Coulomb interaction, shown as a function of cos8 qr at a 



couple of values of the dimensionless parameter, x = rc/r, where rc = 2fiZiZ2e jq . The structure around cos9 qr = —1 is 
due to the deflection of classical trajectories directed towards the center of the repulsive Coulomb interaction. As the relative 
momentum within the pair increases, at fixed r, the range of backward angles where the kernel is significant shrinks. 



B. Kernels for Repulsive Coulomb Interaction 

The Coulomb final-state interaction also provides means for determining S(r) from lZ(q) [IjJ In the Coulomb 
case, there is no straightforward way to invert the relation to the source function as in the case of pure interference. 
To understand how S impacts 1Z for the Coulomb interaction, the structure of the kernel must be considered in detail. 

In the classical limit of a repulsive Coulomb interaction, the kernel depends solely on 9 qr and on the ratio of the 
pair relative-energy to the pair- Coulomb-energy at emission: 

r qZ r 

where \x is the reduced mass and rc is the radius of the Coulomb barrier in the pair interaction. Quantal effects 
become important for qr/h of the order of unity. In the quantal situation, the kernel also depends on the dimensionless 
variable qr/h. The Gamow-factor parameter is then given by the parameter product 77 = xqr/2 = Z\ Z^ e 2 /v and a 
further parameter multiplication yields x (qr) 2 = 2r/ao, where ao = 1/Z\ Z2 e 2 /i is the Bohr radius. 

In the classical limit, the squared wavefunction represents the ratio of initial and final spatial densities during 
motion to the detectors for particles of a given initial momentum, or, due to the conservation of phase-space density, 
the inverse ratio of the momentum densities at a fixed initial separation: 

\A-)i a \|2 d3r °° d 3 9o 1 + cos 6 qr - X 

\4> [ >{q,r, cos 9 qr )\ -> = = === o e(l + cos6> gr -2a;). (17) 



r d 3 q ^/(l + cose qr - xf 



Here, q represents the relative momentum of the particles at separation r, while Too represents a remote separation 
of the particles at some instant when the relative momentum has approached q. The r.h.s. result in (JTTJ) follows from 
Coulomb-trajectory considerations [29(. The implications of this result simplify |25j for £ = and in the high-energy 
(x <C 1) and low-energy (x ^> 1) limits. Thus, for £ = 0, the partial kernel is 

M^) = 4^-i = *°-i = -P=-i. (18) 

<T dg q Vl-x 

The three-dimensional kernel K, is illustrated in Fig. [T] as a function of cos8 qr , for two x values. In the high 
energy limit of x -C 1, the wavefunction squared \(j>\ 2 is close to unity for most values of cos9 qr . In that limit, the 
trajectories emerging from r are practically straight lines. The exception are the trajectories aiming towards the 
center of interaction, that cannot penetrate the region of radius rc- The trajectory directions with cos9 qr < — 1 + 2x 
get shadowed. (The factor of 2 in front of x stems from the fact that trajectories get deflected both during the 
motion towards and away from the barrier.) The deflected trajectories pile up primarily outside of the shadowed 
region. The pile-up produces an integrable singularity in the kernel as shown by the peaks in Fig. [1] At low £, 
the Legendre polynomials have a limited angular resolution. Given that limited resolution, for the purpose of the 
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angular integration in the angular structure in \4>\ 2 , in vicinity of the backward direction cos6 qr = — 1, may be 
approximated by a 5-function multiplied by the integral over the structure: 

^^{q.r.coaOqr)^ » 1 - - 5(1 + cos6 qr ) , for x ■ max (1, £) -C 1 . (19) 

Since Pg(— 1) = (— 1)^, the high-energy kernel can be then approximated as 

Kt(q,r) ~ {-lY +1 x/2, for x < l/max(M) and qr > max{l,£) . (20) 

For more implications of the classical limit (fT7|) for kernels /C^, see ref. 25]. The classical limit of the Coulomb 
interaction is of interest in the context of emission of intermediate-mass fragments (IMF) from low-energy central 
nuclear reactions [3(J. The practical condition of applicability for the classical approach is that QR 3> max(l,£), 
where Q is the characteristic momentum which yields rc ~ R- With Z ~ v4/2, this produces, at low £, the condition 
(e 2 A 3 raw i?) 3 / 2 ^ 2. Given a typical value of R ~ 5 fm for heavy ion reactions, this amounts to A 3 ! 2 ^ 5 and the 
classical limit is met at the fragment mass of A > 6. IMF correlations are commonly expressed in terms of reduced 
velocity v rc d = qj [i(Z\ + Z2) 1 / 2 . With Z ~ A/2, rc ~ 2e 2 /niNVrcd and i? ~ 5fm, the high-energy limit of a: <C 1 for 
the kernel is reached at w ro d <^ 0.04 c. 

At small qr in an interacting system, quantum effects become important as phase-space derealization produces 
diffraction and tunneling. The diffraction affects the angle of an emerging particle at the level of 1/qr, smearing out 
the angular structure in the kernel. Given the resolution of Legendre polynomials changing with £, the diffraction 
effects become consequently important for qr < max (1, 1), with the kernel Kg. then being reduced on average. The 
tunneling into the region of x > 1 becomes significant down to r = 0, for 2rj = xqr = qrc 1- For singly charged 
particles, such as identical pions, this corresponds to v > 0.02 c. Unlike the classical expression (fTT)) . the Coulomb 
wavefunction squared 

\4> { -\q^ cos 9 qr )\ 2 = G(r))\M{i V ,l,-iqr(l + cos 9 qr )\ 2 (21) 
is analytic around r = in the argument r(l + cos8 qr ) = r + rq. In the quantal expression, 

G ^ = (22) 

is the Gamow factor. The terms in the Taylor series for the squared wavefunction (|2"Tj) behave as r l (1 + cos 6 qr ) 1 . Since 
the highest power of cos 9 qr for an £'th term is £, the quantal Coulomb kernel Kg from ([6]) behaves around r = as 
r , similarly to the identity- interference kernel (|12p . The power behavior, associated with the quantal derealization, 
holds up to r ^ min (ao, q^ 1 )- The behavior extends past the classical Coulomb barrier when rc ^ min (ao, q^ 1 ), 
leading to a further dampening of ICi at £ > 1 in the barrier region. 

As an example, in Fig.[2j we consider K-i kernels, for the pK + interaction. The kernels are shown for £ = (0 — 3), at 
relative momenta of q = 15 (left panels) and 75 MeV/c (right panels), as a function of distance r. Regions of higher 
and lower relative momenta within a correlation function lZ(q) are, typically, used to determine, respectively, about 
the short and long-rage features of the source S(r). The solid lines, dashed lines and symbols, represent the classical 
Coulomb kernels, quantal Coulomb kernels and quantal kernels with inclusion of strong interactions. 

We first address the Coulomb kernels at 75 MeV/c. At large r in the right panels of Fig. [21 the classical Coulomb 
kernels ICi are, up to a sign, independent of £ and fall off as 1/r, consistently with Eq. (I20p . Differences emerge at 
low r. At £ > 1 the different classical kernels JCg switch sign (£—1), number of times, cf. [25| . With the distance scale 
of rc being, though, a fraction of a Fermi, rc = 0.17 fm at q = 75MeV/c, the structures associated with the sign 
switching are all pushed to very small r. Those sign-switching structures cannot be seen in the figure as the kernels at 
r < 1.5 fm are suppressed for the sake of clarity. For the quantum kernels, given q^ 1 = 2.6 fm, those kernels are seen 
to fluctuate with r around the classical kernels at r ^ max(l, £) q -1 . At lower r and £ > 1 quantal effects dampen 
the kernels, and in particular, at very low r the classical sign-switching in the kernels gets completely erased with the 
kernels behaving as r l . With inclusion of the effects of strong interactions in the kernels, additional oscillations are 
observed with r within the figure, as the asymptotic behavior of the partial wavefunctions changes for low angular 
momenta, producing wave-function distortions at moderate distances. We will focus on strong interactions in the 
next subsection. 

At q = 15MeV/c, the classical Coulomb radius rises to rc — 4.1 fm, on one hand pushing the asymptotic high- 
energy region out of the r-range shown in Fig. [21 and, on the other hand, making some of the low-r details of the 
classical Coulomb kernels visible within the shown range. Also, since g _1 = 13.1 fm then, much of the region for the 
left panels of Fig. [2] is dominated by quantal effects. The kernels with inclusion of strong interactions get to be very 
close to the Coulomb kernels, since the phase shifts, characterizing changes in the partial-wave asymptotic behavior, 
are close to zero at low q. 
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FIG. 2: (Color online) Kernels tCt, multiplied by (— 1) , for pK + interactions at q = 15MeV/c (left panels) and q — 75MeV/c 
(right panels), shown as a function of 1.5 fm < r < 30 fm, for £ < 3. Classical and quantal Coulomb kernels and quantal kernels 
with inclusion of strong interactions are represented by solid lines, dashed lines and symbols, respectively. The vertical scale 
may change from a panel to panel to emphasize details. Kernels are suppressed at the lowest r for the sake of clarity. 



The quantal suppression of the kernels at qr < £, progressing with £, makes it difficult to learn about the directional 
characteristics of a source from the particle correlations at low q < = 40MeV/c for R = 5fm. Small kernels 
produce small contributions to the correlation which are difficult to detect experimentally. At the other end, at 
high qr, while the kernels become independent of £, they also become small, falling off as l/(q 2 r), cf. (f20|) . As a 
consequence, an optimal bracket of q for accessing the directional characteristic of sources with typical size emerges, 
such as q ~ (30 — 100) MeV/c for the pK + system. 

We have seen the importance of the analyticity of \4>\ 2 in r around r = 0, which implied the r e suppression of the ICg 
(£ > 1) kernels at low r. The analyticity in q is of interest in the context of the analysis of correlation functions as a 
function of q. In the case of pure identity interference, the wavefunction squared (1 1 1 [) is analytic in q which makes the 
kernels from ([3]) analytic in q. The resulting correlation coefficients TZe are also analytic and behave as q e at small q. 
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In contrast, because of the general dependence on the parameter 



V = , (23) 



inversely proportional to q, and because of the appearance of the Gamow factor (|22|) . the square of the Coulomb 
wavefunction (f2"Tj) lacks analyticity in q around q = 0. Still, it may be tempting to factor out the Gamow factor, as 
has been done in the analysis of pion correlations, arriving, with \(j>\ 2 /G = \M\ 2 , at a Gamow-corrected kernel. The 
corrected kernel remains finite at a given r and 9 rq , when q — > 0, but retains a dependence in this limit on cos9 qr for 
any finite value of r. Thus, \M\ reaches a different value depending on the side from which q = is reached and lacks 
therefore analyticity in q. Correspondingly, even the Gamow-corrected correlation functions from ([3]) lack analyticity 
in q. The Gamow-corrected kernels ICi tend to finite values as q — > 0, rather than behaving as q' . 

As to the high-g limit of the Coulomb kernels, while the resolving power becomes independent of £, the accessed 
information on the source represents one moment only per correlation coefficient, as 

q— »oo q z \r / 

c g— >oo g z \T / 

where 

(F) = /" d 3 rS(r)F(r). (26) 

Since the classical approximation is accurate at large qr, this relation is unchanged by quantum considerations. 
Working within the classical approximation, one might be tempted to continue the tq jt kernel expansion, producing 
expansion of the correlation in powers of 1/q 2 - The next expansion term could give access to the moment of the 
source function with a factor of 1/r in place of 1/r. However, beyond the leading term, the expansion gets altered 
by quantum effects making the expectation invalid for lighter particle pairs such as pK + . 



C. Kernels for Short- Range Interactions 



Strong interactions between two particles also provide leverage for extracting source information from measured 
correlations and are particularly effective when the strong- interaction cross-sections are large [29(. In analysis of the 
effects of the interaction, the two-particle scattering wavefunction may be conveniently decomposed into angular- 
momentum eigenstates, here with any intrinsic spins suppressed: 



%r,cos V) = E a ^ P ^( cos 9 v) e i?rcose - +/* (tt - 9 qr ) — . (27) 



r — >oo 



The coefficients ae in front of the radial wavefunctions Ri are adjusted to make the large-r outgoing wave-contributions 
to <^( - ) identical to those for a wave in the absence of short-range interactions. The explicit asymptotic form of the 
wavefunction, on the r.h.s., is shown for the case with no Coulomb interactions. The amplitude / = is the 

scattering amplitude associated with the more standard scattering wavefunction (f>^ where conditions are imposed 
onto incoming wave-contributions in the asymptotic region. The amplitude factors in the two wavefunctions are 
related with f^(9) = [/<+>(* -6)]*. 
The three-dimensional kernel from (|27|l is 

K(q, r, cos 9 qr ) = |» - l = _ ^ 0/i a \ R tl (r) R l2 (r) P tl (cos 9 qr ) P (2 (cos 9 qr ) - 1 (28) 

T tit* 

Since most shape information is contained at intermediate to large r, the large r limit is especially insightful and can 
be found from Eq. Ij2"7]). 

IC(q,r,cos9 qr ) w \ ^-(ir - 9 qr ) - Im — — — sin[gr(l + cos 9 qr )] 
r—>oa r ail r 

+ Re/(7r-M cos[gr(1 + CQS Qqr)] (29) 
^ ^^-M- £^(1 +cosM- (30) 
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FIG. 3: In the asymptotic region of r — > oo, for short-range interactions, the kernel K represents changes to the flux of particles 
emerging at relative momentum q. produced by scattering. Here, the particles were initially separated by r. 



In arriving at (pZ9"|) and ([50]) . we used an expression for the scattering cross section, |/| 2 = da/dfl, and the optical 
theorem, Im/(0) = qa/Air. Further, in obtaining the expression (|30[) we have exploited the fact that for qr — > oo, 
under an r-integration, the respective products in (|29[) approach the limits: qrsin(qrx) — > S(x) and qr cos(qrx) — ► 0. 
The result (|30[) has a simple geometric interpretation. The first term represents the probability that a particle aimed 
directly at the scatterer reflects into the direction of q, while the negative term represents the shadowing of a spherical 
angle, a/Airr 2 , around the backward direction. 

The partial kernels fCi, respectively from Eq. (|2"5)) and ([50]) . are 

£i(q,r) = jXifo 0* oj a ti a *iz R ti{ r ) R iA'r)-ho (31) 



(-1) 



(32) 



where we have made use of the Wigner 3j-symbol. In the asymptotic region, all £ > 1 kernels decrease as 1/r 2 and 
are proportional to the scattering cross section. Larger cross sections, such as associated with resonances, will lead to 
larger kernels at a given q. Otherwise, the kernels are sensitive to the angular dependence of the cross sections. For 
q in the range for which correlations are typically studied, cross sections are usually fairly isotropic. Thus, the first 
term in the brackets on the r.h.s. of 1|32[) decreases quickly as £ increases, yielding 



JCe(q, r) sa — - — - — , for large r and £. (33) 
iirr 

Apart from the sign, this kernel is independent of £, similarly to the Coulomb kernel (|2"U|) . In analogy to the classical 
Coulomb results for individual £, the asymptotic results (|32j) require qr 3> £ for their validity. An additional condition 
for validity is that r/d^> £ where d is the interaction range. 

Regarding the limiting short-range (|33f> and Coulomb (|20p /Q-results, their independence from £ results from the 
kernels being exclusively influenced by the shadowing of particles by each other. The effect of shadowing decreases 
with distance between the particles, but more slowly for the Coulomb than for purely short-range interaction. In the 
case of short-range interactions, the kernel falls off as <r/r 2 , while for the Coulomb forces the kernels fall off as rc/r, 
where rc = 1liZ\Z%^ / q 2 is the radius of the Coulomb barrier. Since rc explicitly decreases with q, Coulomb kernels 
decrease monotonically with increasing q, whereas for the strong interaction, the q dependence on the kernel tends to 
follow the energy dependence of the cross section. Strong interactions can thus provide resolving power at large q, if 
the cross sections are large, e.g. at resonance energies. 

When the two terms in the asymptotic form of the kernel (|30p are integrated over full spherical angle, those terms 
exactly cancel, reflecting the fact that, for a short-range interaction, the particles which are shadowed for cos9 qr = — 1 
reappear at other angles at the same magnitude of relative momentum. Consequently, the £ — coefficient for the 
1/r 2 asymptotic fall-off in (|3li|) vanishes, implying that the /Co kernel from (|3"Tj) vanishes faster than 1/r 2 as r — * oo. 
These results may be contrasted with what is found for the Coulomb interactions. The large-r kernel Coulomb kernel 
from (flUl) yields a finite result when integrated over the spherical angle. In direct consequence, the £ = kernel is 
finite in the asymptotic region, Eqs. (JTHJ) and (|2"U|). The difference, between the short-range and Coulomb cases, is due 
to the fact that a Coulomb interaction does not just change the orientation of the relative momentum but also makes 
the momentum magnitude different from the asymptotic value at any finite distance from the interaction center. 
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At low relative momenta q, s-wave scattering is likely to dominate. An important result, visible already in the 
asymptotic limit, Eqs. (|32p and is that the s-wave scattering contributes to the kernels at all £. At the level 

of Eq. (f3"Tj) . it is seen as associated with the interference between the waves. At resonance, the radial wavefunctions 
Ri will be particularly large in the near zone giving rise to enhanced Kg. there, including £ = 0. At low q and finite 
£, the effects of the centrifugal barrier can be important. When q is so low that the barrier is encountered outside 
of the strong-interaction range, qr < I + i at r > d, the radial wavefunctions are approximately given in terms of 
the spherical Bessel functions, Rt/r oc je(qr) oc (qr) . With i?o approaching a constant value as a function of both 
q and r in the region d < r <C 1/g, at q — > 0, given the triangle inequality for the £- values in the 3j symbol in (132p . 
the strong-interaction kernels will behave as Kg(q,r) oc (qr) l (l + b/r) for £ > 1 in the region. The kernels will be 
suppressed then at low q and low r as well for larger £. A short-range repulsion may enhance suppression at lowest 
r < d beyond that produced by the centrifugal barrier. 

Outside of the interaction range, at r > d, the wavefunctions Rg are fully determined by the respective strong- 
interaction phase-shifts, both in the absence and presence of Coulomb interactions. We exploit this [291 ] in obtaining 
the kernels for strong interactions, e.g., the kernels in Fig. Energy derivatives additionally constrain the integrals of 
the wavefunctions Rg squared over the region r < d [29L l3l| . This facilitates an extrapolation of the kernels down to 
r = 0, adequate for source resolutions that can be practically achieved, without the need to resort to potential models 
for the strong interaction. Still, with the strong-interaction kernels becoming model-dependent at short distances, 
we refrain from displaying them at r < 1.5 fm. Needless to say that, given the quark substructure of hadrons, the 
description in terms of wavefunction for relative motion becomes questionable at the shortest distances r ^ 0.5 fm. 
For IMFs the same becomes true at even larger distances, due to the finite size of the nuclei. 



III. SPHERICAL CARTESIAN HARMONICS 
A. Harmonics within Spherical Coordinates 



We start out with a review of basic concepts concerning surface spherical harmonics. A harmonic function is one 
which satisfies the Laplace equation 

V 2 F(r) = . (34) 

When using spherical coordinates, it is convenient to express harmonic functions, which are regular at the origin, in 
terms of tesseral harmonics Yg m , 

F(r ) = ]T F; m t 1 Y lm (O) , (35) 

Im 

where the coefficients Fi m take on arbitrary values. In the summation over t, the individual terms in (|35p . of the form 

(r) = V^r e Y, Fim Y tm (Sl) , (36) 

m 

are harmonic functions of degree £. The latter stand for homogenous functions of degree £ in x, y and z, that are 
harmonic. A surface spherical of degree £ is a harmonic function of degree I, taken on a unit sphere, r = 1: 

FW(fi) = VtoJ2*tm Ytmffl . (37) 

m 

Given the completeness relation for the tesseral spherical harmonics, 

s(n - n') = -L ]T(2£ + 1) p e (n • n') = £ r; m (fi') r, m (o) , (38) 

l im 



any function of spherical angle, and in particular 1Z and S in Eqs. (|8|)-(|10|). can be expanded in the tesseral harmonics 
upon using (|38p with the identity 

G(fi) = / dfi' S(fl - Q') G(n') . (39) 



Equation (|35p for the harmonic functions represents an example of the expansion, where the expansion coefficients 
have a specific radial dependence. Within every rank £ of an expansion, the expanded functions are generally described 



11 



in terms of (21 + 1) complex coefficients Fg m . However, for real expanded functions, the coefficients F^^ m= o must be 
real and, otherwise, the coefficients must satisfy F^_ m = (— l) m F/ m . Thus, for each £, the expansion is described 
in terms (21 + 1) independent real numbers. Under rotations of the coordinate system, for a given I, the tesseral 
functions transform like the components of a spherical tensor of rank I. For a function that is independent of the 
coordinate choice, the expansion coefficients transform further in this fashion; the sums over m in (|35|) . (|36[) and ([5]) 
represent scalar products of the spherical tensors. In fact, the r.h.s. of Eq. (f38|) also represents a superposition of 
the scalar products and, in this, represents a covariant generalization of the middle result in (|38p . which is seen as a 
superposition of the scalar products of the tensors in the case when the coordinate z-axis is oriented either along n 
or n'. 



B. Cartesian Harmonics 

Cartesian harmonics allow functions of the spherical angle to be expressed as superpositions of the scalar products 
of cartesian tensors. For real functions, the coefficients of expansion in cartesian harmonics are real, which is not 
necessarily the case when expanding functions in tesseral harmonics. The expansion coefficients for cartesian harmonics 
are generally easier to interpret than those for tesseral harmonics. 

Let us consider a homogenous function of degree £. When using cartesian coordinates, a homogenous function may 
be represented as 

F^\r)= ^ ' Faia 2 . ..a e r ai r a 2 ■ ■ ■ r a, > (40) 

Ql Q2 ■ ■ 

where a = x,y,z are cartesian indices. The following remarks can be made regarding this representation in terms 
of cartesian coordinates. Since the product of cartesian coordinates on the r.h.s. of (|4"0")) is symmetric under the 
interchange of its terms, attention may be limited to the coefficients of expansion F ai0l2 „_ ai which are symmetric 
under the interchange of indices. Since the product of t coordinates transforms further as a cartesian tensor of rank 
£, the coefficients F ai a2 ... af , must transform as a cartesian tensor of rank £ for a function F^ that is independent of 
the choice of directions for coordinate axes. Given the symmetry of the tensorial coefficients, the r.h.s. of (|4"0"]) can 
be conveniently rewritten by grouping terms with different numbers, £ x , £ y and t z , respectively, of the x, y and z 
coefficients, obtaining the representation such as in ([5]), 

F(t) 00 = E 777777 V" y iv **' ^ T ^ ( 41 ) 

i 

where £ = £ x + £ y + l z . The last expression in (|41[) is a schematic representation for the convolution of the cartesian 
tensors. 

On substituting the representation (|40[) into the Laplace equation (|34|). we find that the homogenous function is 
harmonic if and only if the cartesian tensorial coefficients are traceless [12|, i-e., 

^ ] ^ail..tt t -2 a a ~ @ I (42) 

a 

or 

Fi x +2,e v ,e z + Ft x ,ty+%,i z + Ft x ,e y ,t x +2 = ; (43) 

where £ x + £ y + £ z + 2 = £. The components for a symmetric tensor, labelled with the values of (£ x ,£ y ,£ z ), are 
illustrated with a triangle in Fig. Q] When moving along lines parallel to the different sides of the triangle, either 
the value of £ x , £ y or £ z stays constant. The tracelessness condition (1431) relates components at the corners of an 
equilateral triangle two £-units on the side in the figure. 

In the appendix we show that there exists a real projection operator V, which operates in the space of cartesian 
tensors of rank £ and projects, out of any tensor, the part which is symmetric and traceless. As an operator in the 
space of rank £ tensors, V itself is a rank 2£ tensor, £ times covariant and £ times contravariant. As a projection 
operator, V satisfies V 2 — V and "P T = V . In terms of the operator V, for any tensor T of rank £, the function given 
by 

(VT) T r e = T T V T r e =T T Vr e (44) 
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FIG. 4: Triangular diagram representing the components of a symmetric tensor, in terms of dots labelled with the component 
[l x ,£ y ,£ z ] values, for t = £ x + £ y + £ z = 4, after [261 ] . When moving along the lines parallel to the respective sides of the triangle, 
either l x , £ y or £ z stays constant. The tracelessness condition (|43[) relates tensor components at the corners of an equilateral 
triangle two £-units on the side, within the diagram, such as those marked by the squares. As the set of linearly independent 
components of a traceless symmetric tensor, the £ z — and i z = 1 components may be chosen, marked with the thick grey 
lines in the figure. 



is harmonic, where we follow the schematic notation. Since, the function (|44|) is harmonic for any tensor T, each of 
the components of Vr l on the r.h.s. must be harmonic. On the sphere of r = 1, those components define the cartesian 
surface spherical harmonics 

•^l]a2...a,(^) = E ^1 a 2 —a e :a' 1 a' 2 ...a' e na 'i n °2 ' " ' ""J ' 0^) 
a[ a' 2 ...a' t 

where n — (sin 8 cos <fi. sin 6 sin 0, cos 9), which can be further also written as 

Mty = E 77777777 V M> n v ^ ■ ( 46 ) 

V u 

Cartesian components of the direction vector n are combinations of the tesseral rank-1 harmonics Y\ m - The product 
of £ components, as e.g. demonstrated by the quantal rules of angular-momentum superposition, is a superposition of 
tesseral harmonics of rank I and lower, of the same evenness as I: 



bmcc are harmonic, can consist of the tesseral rank-£ harmonics only. As will be demonstrated in the 

appendix, the operator V consists out of a symmetrization operator, i.e. the unit tensor in the space of symmetric 
cartesian tensors, and of operator terms which involve taking traces of one or more pairs of the cartesian indices. 
In acting on the the cartesian product n e , the symmetrization operator leaves the product n l intact while the other 
operator terms in V produce symmetrized terms of the form n e ~ 2k {S) k , where k — 1, . . . ,1/2. The latter terms 
represent superpositions of tesseral harmonics up to the rank of t — 2, of the same evenness as t. The implication is 



E 

I'm' 
£'<£ 



I'm' ^ 



Ygi m r(£i) 
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TABLE I: Cartesian harmonics for I < 4. Other harmonics can be found by either permuting the indices, i.e. A xy x = Axxy, or 
by swapping indices on both sides of an equality, e.g. x *-* y. Thus, given Axxy = n 2 n y ~ 0-/5)riy, swapping the indices y <-> z 
yields Axxz =n%n z - (l/5)n 2 . 





Tlx 


/i( 3 ) _ 

<s\.xyz — T^x f^y l^z 


A {2) - n 2 


-1/3 


A { xLx =n%- (6/7)n 2 x + 3/35 


A ( xy ] = n 


x ^y 


■Axxxy — W>x ^y (3/T)Tlaj Tly 


A {3) -n 3 

-t^-XXX ,L x 


{3/5)n x 


Ai%y =n 2 x n 2 y - (l/7)nl - (l/7)n 2 y + 1/35 


A {3) -n 2 n 

<r\-xxy — lb x lb y 


- (1/5)% 


■A-xxyz — n>x ^y ^2 \\i'^j'^ / y ^ 2 



that the projection operator V removes from n e a portion represented by the superposition of the lower rank tesseral 
harmonics (in parenthesis on the r.h.s. of (|47[) ) while retaining the highest rank £ portion intact. Given the general 
form of A^\ from application of V onto n e , the cartesian harmonics may be constructed recursively [25|], based on 
the tracelessness of the A^ tensors within each rank, starting with _4^ =0 ) = 1 : 

A W = iVn A {1 ' 1] 

i=l 

~ I (11 — 1) ^ 1 ^ ' ^jiij n a A a a 1 ---a i - 1 a i+1 ---a j - 1 a j+1 ---ai ■ (48) 
^ ' l<i<j<£ a 

Expressions for cartesian harmonics of lowest rank are given in Tableland the recursion with (|48p . otherwise, produces 
the series ([7]). One of the features of the cartesian harmonics is that the three cartesian axes are treated equally, 
unlike in the case of tesseral harmonics. Within a given rank £, not all cartesian harmonics are independent. A 
symmetric tensor has (£ + 1)(£ + 2)/2 different components, cf. Fig. 5] The tracelessness condition (I43[) reduces the 
number of linearly independent components, including the number of independent A^ functions, to (2£ + 1), equal 
to the number of Yg m functions. As a linearly independent set, one might choose the components with £ z = and 
£ z — 1, see the figure. 

In the case of ot\ = a.2 = ■ ■ • = ote = z in (|47p . both the cartesian component product n l and A^> are invariant 
with respect to rotations around the z axis. For invariance, both functions must be a superposition of m! — tesseral 
harmonics only. It follows then that A^kl.z must be proportional to Yg Q or Pp. Since the action of the projection operator 
V leaves the highest rank terms unchanged, the coefficient for Pi in A^ must be the same as in n e . Accordingly, the 
respective cartesian harmonics are given by 

4L.,(")= pj^ ^(coBfl). (49) 

Starting with the first equality in (|38p . we are now in position to produce a completeness relation in terms of the 
cartesian harmonics. Thus, if n' is directed along the z-axis in (j3"8|) , we can express, given (pf9"|) , the Legendre 
polynomial in the middle term in (f3"8")) as ((2£ — l)ll/£\)(n' e ) T A. When rotating n' and n, to move n' away from 
the alignment with the z-axis, the scalar product of the tensors transforms covariantly, allowing to obtain different 
general representations for the 5-function in spherical angle in terms of the cartesian harmonics 

5(Sl-tf) = — ^(2^ + 1) Pe(n-n') 



Ql...Qf 



i E (2i+i)u E A^w^um 



Ql...Qf 



I 

= h^J^MV)<*<*n{». (50) 

x V z 

In obtaining the third equality in (|50[) and later, we use the projection properties of the operator V: 

(n' e ) T A= (n' e ) T VA= (Vn ,e ) T A = (A') T A. (51) 
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Use of (|50|) in the identity (|39[) allows for functions of spherical angle to be expressed as a sum of cartesian harmonics, 
such as in ([5]). Independent of the orientation of the coordinate axes, the coefficients of expansion in cartesian 
harmonics, such as in ([9|), transform as components of a cartesian tensor under rotations. As the tensor is traceless, 
only {2£ + 1) coefficients are independent within each rank I. 



C. Operations with Cartesian Harmonics 



When given some specific expansion coefficients in terms of cartesian harmonics, it may be of interest to calculate 
the associated cartesian moments. To facilitate that, it is useful to express products of cartesian components in terms 
of cartesian harmonics, which can be done [25| through recursion starting from Eq. obtaining 



E 



0<mi<£i/2 
X 



{2l-Am + l)\\ 
2 m (2£-2m+l)U 

£ x ! 



f 1 
' v ■ 



(4-2m x )!m x ! {£ y -2m v )\m y \ {£ z -2m z )\m z \ e ~ 2 ™ 



(52) 



The rank of A within each cartesian component on the r.h.s. is equal to or lower than the power of the corresponding 
cartesian component on the l.h.s. From ([52")) and ^j, we find 



4ir£l ^ {2£-Am+l)\\ 
(21+ 1)!! ^ 2 m (2£-2m + l)\\ 

rh 

0<mi<^/2 



(53) 



Ly. 



4! 



/ drr^ 2 ^„ M (r), 
Jo 



(£ x -2m x )lm x \ (£ y - 2m y )\ m y \ {£ z -2m z )\m z \ J 

which allows the moments to be obtained from cartesian harmonic coefficients. 

As combinations of tesseral harmonics of different rank £, the cartesian harmonics of different rank are orthogonal 
to one another. Due to different behavior under inversion, see Eq. ([7]), two cartesian harmonics are further orthogonal 
if any pair of their li subscripts differs by an odd number. Otherwise, the scalar product of harmonics is given by (2r| 



W) 



47T 



{21 +1)U(2£- 1)!! 



E 



to! (2l-2m-l)!! 



n 

i=x,y 1 z 



max(0,(£--£i)/2)<mi 
iTii<min (e' i -£ i /2,e i /2) 



((£' i -£ i )/2 + m i )\(£ i -2m i )\m i \ ' 



(54) 



The anticipated symmetry between £ and £' on the r.h.s. of (|54|) may be seen by eliminating TOj in favor of m! i — 
mi + (4-4)/2. 

Cartesian expansion coefficients may be converted into tesseral expansion coefficients and vice versa. In terms of 
cartesian harmonics, the m > tesseral harmonics are given by [26| 



Y lm (Q) = {-l) rn {21 - \)\\ 



2£ + l 



_4tt {£ + m)\ {£ - m)\ 
The cartesian harmonics may be, on the other hand, expressed as 



1/2 m 
k=0 



■A.\u.: 



(m-fc)!fc! [k 



-k,e- 



(55) 



(56) 



where 26], for m > 0, 

of = (-l) m (2^-1)!! 



2£+ 1 



4tt {£ + m)\{£-m)\ 



1/2 m 



E-' 

fe=0 



i — k 



(to — fc)! fc! 



(£\[k,m-k,£-m\) 



(57) 
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TABLE II: Transformation equations between cartesian and tesseral harmonics for I < 4. 
Yio = (1/2)^37^ Ai 1} 

Yi±i = -(i/2)y/3/2^A^ T (1/2) 73/277 4 X) 
Y 2 o = (3/4) 75/ 77 4 2) 

Y 2 ±i = -(V2)VW2^4 2 i T (V2) y /TEfa A 2 } 

Y 2 ±2 = -(1/4) yiHT^ A$ ± (i/2)y/TK/2^A^ + (1/4)^/1572^ 42 

Y i0 = (5/4) 7^777 4 3) * 

y 3 ±1 = - (5t/8) y/217^4 3 ^ T (5/8) 7 2T7774 3 ) z 

y 3 ±2 = -(1/4)^/105/27^^ ± (i/2)V 105 / 27r ^ + (l/4)x/105/27r^ 3 2 z 

y 3±3 = (i/8)VW^4ra ± (3/8)^/3577?^^ - (3i/8) 735777 47i, =f {l/8)y/35/nA£L 

y 40 = (105/16) jlf^A^L 

y 4 ±i = -(21^/8)^57^4^ T (21/8) V / 57ir-4 4 i« 

y 4 ±2 = -(21/8)^5727? 4t« ± (2H/4)V572^4 4 !/ ) ^ + (21/8)^/5727741, 

y 4 ±3 = (3i/8) A v % z ± (9/8) 735777 41L - (9i/8) 7^ 7^ 4« U T (3/8) 735777 41, 

y 4 ±4 = (3/16)V35/27r4tU T (3i/4) V35/27T 41]^ - (9/8^/35/2^4^ ± (3i/4)735/27r 417 y + (3/16)735/2tt .4! 
4 1} = 2v7773yio 

4 X) = 4^2^/3 yi +1^/277/3 y_i 

4 1} = - y /2lF/3Y 11 + ^lTJlY X -^ 

4? = (4 /3) 7777 5 y 20 

4* = VWiS^i +i/27r/15y 2 _i 

4? = -\/2777T5 y 22 - (2/3)777/5^0 - 7Wi5y2- 2 

4 2) = -TwT^i + 72 *715* 2-1 
4? = - V 2 ^/i5y 22 + i7 27 i-/i5y 2 - 2 

4 2 2 = 727715 y 2 2 - (2/3) 77775 y 20 + ^JTEy 2 -2 
a z % = (4/5) 77777 y 30 

4« = (4^/ 5)777/2 1 y 3 i + (4i/5)77772Ty-i 

4 3) * = -727r/i05y 32 - (2/5)77777^0 - ^/ToEy^ 2 

■Ayyy = -i 7*735*3 3 - 7/5) 7^77 y 3 1 - (i/5) 737777 y 3 -i -i777735 y 3 - 3 

-(4/5) 777/21 Y31 + (4/5)77772Ty 3 -i 
4 3 Jz = -i727r/105y i2 + i727r/105y 3 _ 2 

4^ = 777735^ 3 + (1/5)^721^1 - (1/ 5) 7^72 1 y 3 _i - 777735 y 3 _ 3 
4 3) z = 727r/i05y 32 - (2/5) 77777 y 30 + 72777105 y 3 -2 

Ai% = i 7^735 Y 3:i - (i/5) 77772T y 3 i - (i/5)7*721*3-i +17*735*3-3 

4 3 i = -777735 y 33 + (1/5) 73777 y 3 i - (i/5)737777y3-i + 777/35 y 3 _ 3 

A z % z = (16/lO5)0Fy 4O 

A ( yL = (4i/21) 77775*41 + (4i/21) 777775*4-1 

A y % z = -(2/21)727775 y 4 2 - (8/105)7^40 - (2/21)727775 y 4 -2 

A y % z = -7/3) 777/35 y 43 - 7/7) 7775 y 4i - 7/7) 777/5 y 4 _i - 7/3) 77735 y 4 _ 3 

■Avvvv = (l/3)727r/35y 44 + (2/21)72tt/5 Y 42 + (2/35)7*^4 + (2/21)72tt/5 Y 4 - 2 + (1/3)72tt/35 Ya-a 

41L = -(4/21) 77775 y 4i + (4/21)77775 y 4 _i 
A%z = -(21/21)727775 y 4 2 + (2i/2i) 72775 y 4 _ 2 

4 4 yU = (1/3) 777735^ 3 + (i/2i)77775y 4 i - (1/21)77775 y 4 _i - (1/3)777 /35 y 4 _ 3 

Ai% v = 7/3)7277735 y 44 + 7/21) 7177/5 y 4 2 - 7/21)72*75*4-2 - 7/3)72*735*4-4 

4tL = (2/21)72775 y 4 2 - (8/io5)0Fy 4O + (2/21) 7^775 y 4 _ 2 

Ai%z = 7/3) 07 /35 y 43 - 7/21)7*75 y 4 i - 7/21) 77 75 y 4 _ i + 7/3)7*735 *4- 3 

4x5,5, = -(1/3)727735 y 44 + (2/io5)7¥y 40 - (1/3)727735 y 4 _ 4 

4tL = -(1/3) 7777^ 43 + (i/7)v775y 4 i - (1/7)77775 y 4 -i + (1/3)77773 5 y 4 _ 3 

4 4 *U = -7/3) 72*735 y 4 4 + 7/21)721775*42 - (i/21) 72775 y 4 _ 2 + 7/3) 7277735 y 4 _ 4 

4tL = (1/3) 727735 y 44 - (2/21)727775 y 4 2 + (2/35)0ry 4O - (2/21) 72775 y 4 _ 2 + (1/3)7277/35 y 4 _ 4 
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the overlap integral (-|-) is given in (|54[) and a-" 1 = (— l)™ 1 (a™)*- Formulas for conversion between tesseral and 
cartesian harmonics for I < 4 are given in Tablefe] 

Many circumstances require multiplication or division of the functions of spherical angle. For example, the empirical 
correlation function C — 1 + 1Z, as a function of spherical angle at a given P and q, follows from an equation of the 
form 

x(n) = c(n)B{n), (58) 

where the l.h.s. represents the measured two-particle yield, X(Q) = ,j d6j \ 3 — , cf. Eq. JT]), while B on the r.h.s. is the 

Q Pa Q Pi) 

product of single-particle yields, 33^, and possibly of two-particle efficiency In order to find C(f2), one could 
bin the the three functions in two dimensions of cos 9 and <p and obtain C within each bin by dividing the value of X 
in the bin by the value of B in the bin. However, if one next wanted to expand C in cartesian harmonics, the result 
could be distorted by residual binning effects which would disappear in the limit of fine binning, but that could be 
only afforded for high statistics. 

An alternative approach is to explicitly solve for the coefficients Cg given Xj and Bp. With this, one could forego 
the binning in cos# and <j). For instance, when assessing X and identifying a particle pair at an angle f2, one could 
increment values in an array for Xg by \{2t + l)V./£\\AACt) as described by Eq. ©. One could then carry out an 
analogous procedure for the mixed pairs, filtering them through two-particle efficiency. In this way, both Xf and Bp 
could be obtained in a straightforward fashion, without any binning in cos 9 and (f>. 

Though the determination of Xp and Bp could be straightforward, the determination of Cg can be more difficult. 
Here, we outline one possible procedure based on the expansion (J5J) of the functions in powers of q. If we introduce 

we can write 

Jr(g) = E^^"W* = E E 'y(?-?)'y@) B t-i> c r>&&&- ( 60 ) 

If the condition of tracelessness for Cf is lifted, the solution to Eq. (|60p ceases to be unique, because q~i are not 
independent due to the constraint q% + q^ + q% = 1. We will first find a solution to (JBDJ) without imposing the 
tracelessness condition, which will provide us with C(f2) from which we can next find the traceless Cp 

The simplest strategy in finding a solution to (|60p is to treat the components q~i as if they were independent. If we 
denote the solution by Cg, we obtain 

d 

0<l[<l.i 

which may be solved for Cf iteratively. The £ = term follows right away, 

Co = ^ , (62) 



B, 







while the subsequent coefficients follow from 



°* = T?--i¥rZr E l$-h^) B t-eC t ,. (63) 
l'<t—\ 

The combination of 

C(Q) = J27$ (64) 

t 

and of © next yields for the traceless coefficients 

c i = I J °W = E 7(* + 2m) (¥+ 2m] C i+2A , (65) 
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with 

(21+1)!! v 1 4! ^ ! 

(£ + £' + 1)!! V 2 m - 2m x )lm x l (£' y - 2m y )\m y \ (£' z 

0<mi<£'/2 
m=(t'-i)/2 

where in the second equality we have made use of the expansion (|52p for the product of direction vector components. 

In practical circumstances, the sums over I must be cut off at some £ max . This is because, on one hand, the 
experimental apparatus and the statistic limit the angular resolution and, on the other, the features of a reaction 
limit the number of £ for which C^may be significant. Exceeding a sensible ! max will result in coefficients dominated 
by noise. A likely situation is that of a suitable value of £ max higher for X, B and than for C . This is because 
the anisotropics associated with measurement efficiencies, particularly at small q, can be much stronger than physical 
anisotropics, especially when the experiment is blind in certain directions. In the latter case, the binning method 
would necessarily fail for small bins, while the method relying on the ^-decomposition from start could succeed when 
using large £ max for X and B, as long as only low £- values were demanded of C. 

Strategies similar to that employed for finding C in Eq. (|58[) may be employed for finding the expansion coefficients 
for an inverse of an angular function or for a product of angular functions. In finding the inverse of B in terms of B^ 
one would replace X(il) in (f58"| by 1, leading to Xg = 8go in Eqs. (|62p and (|63[) . On the other hand, in looking for 
the product of B and C in terms of expansion coefficients, one can first obtain the coefficients 

X /=4l £ lil-^li^B^Ct,, (67) 

that do not meet the condition of tracelessness. Next, one can obtain traceless coefficients X? following such a 
procedure as for C in Eqs. (|65[) - (|66p . Obviously, that specific detracing procedure is a form of application of the 
projection operator V, described in the appendix, to a symmetric tensor. 

We have, thus, demonstrated some practical procedures for dividing, inverting and multiplying angular functions in 
terms of ^-arrays without resorting to binning of the functions in cos 9 or <j) and without performing matrix inversion 
operations. 

IV. PRACTICAL ASPECTS OF SHAPE ANALYSIS 

The expansion of correlations and sources in tesseral or cartesian harmonics brings about practical benefits. The 
three-dimensional dependencies on relative momentum and spatial separation of the functions are replaced by sets 
of one-dimensional dependencies of angular expansion coefficients. The replacing of the three-dimensional by one- 
dimensional dependencies allows, principally, for a more thorough representation of the correlation data and of the 
physical picture of those data's implications. The three-dimensional integral relation between the correlation function 
and source gets replaced by a set of one-dimensional relations between the respective expansion coefficients, simplifying 
the task of finding the source features. The expansion coefficients for the functions, as either spherical or cartesian 
tensors, transform covariantly under rotations. 

Of the two sets of harmonic functions for expansion, the tesseral set represents an orthonormal basis, but leads to 
coefficients which are complex and hard to interpret. On the other hand, the cartesian set is overcomplete, generally 
nonorthogonal within any rank, but yields real coefficients with relatively straightforward interpretations. The trans- 
formation properties, under rotations, are more straightforward for the cartesian than for the tesseral coefficients. We 
have shown, in the preceding section, how to transcribe results between the two harmonic representations. 

A. Shapes in Terms of Cartesian Harmonics 

Along any direction of relative momentum or separation, the correlation or source function are equal to combinations 
of the cartesian coefficients ((8|). Particularly simple results follow along any of the employed cartesian axes; e.g. along 
the positive direction for the x-axis the correlation function is: 



5' ., r (£\i - 2m) , (66) 



TZ(q) = KM (q) + (q) + (q) + TZ^q) + 



(68) 
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Terms for subsequent I refine the information on the dependence of the function on g, along the direction. As we have 
mentioned, the physics, on one hand, and statistics and apparatus, on the other, will limit the number of significant 
terms in the expansion, with the high-^ terms being dominated by noise (24j . 

For identical particles, the correlation function and accessed emission source are invariant under inversion. In such 
a case, the odd-£ terms in the harmonic expansion vanish. That is, generally, not the case for nonidentical particles. 
For instance, if the protons in a reaction are, on the average, emitted earlier than pions, the distribution S of pairs 
in r p — r n will be pushed out in the outward direction, away from 0. The three £ = 1 cartesian moments (|52l) . 

computed from Sa \ yield the magnitude and precise direction of the average displacement of the distribution. At 
individual r, the effects onto S of different emission times and possibly of different emission locations will result in a 
dipolc distortion representable in terms of magnitude and direction: 

Si 1 Hr) = S^(r)ei 1 Hr). (69) 

In the above, the direction vector is e (1) = [sinflWcos^^sinfl^sin^^cos^ 1 )]. Not only the magnitude of the 
angular distortion, but also the characteristic directions of the distortion can depend on the relative distance. If the 
distortion direction is independent of the distance, this direction will be further the same for the average displacement 
and the same for the correlation at any relative momentum. The latter is due to the fact that the kernel ICi is 
universal within a given rank, i.e. the same for different tensor components. The I — 1 distortions are behind the 
effects investigated by Gelderloos et al. [32j and Lednicky et al. [33|, who have compared pair emission probabilities 
for pairs with relative momenta directed along and opposite to the pair total momentum. The advantage of using 
the cartesian harmonics is the ability to determine the distortion directions in correlation, and the implication of 
those distortions for source distortions, without any presumption regarding the distortion direction when no possible 
symmetry arguments can be invoked. 

Quadratic moments of the emission source are associated with the rank i = 2 and t — expansion coefficients, 
cf. Eq. (f53|) . with I = 2 rank giving rise to an anisotropy of the tensor out of the moments, (r ai r a2 ). A triaxial ellipsoid 
could be associated with those moments, with Euler angles describing rotation from the coordinate axes to the axes 
of the ellipsoid, i.e. the axes along which the tensor of moments gets diagonalized. At any value of the corresponding 
relative variable, for either the source or the correlation function, the rank £ — 2 symmetric and traceless tensor of 
expansion coefficients may be diagonalized and described in terms of 5 independent parameters, 2 distortions and 3 
angles: 

c(2) / \ _ c(2) (2) (2) _ f(? (2) g(2), (2) (2) q (2) (2) (2) ( , 

°qi a 2 V I — S l e l«i e lo 2 Wl + D 3 ) e 2a t e 2 a 2 + D 3 e 3«, e 3o 2 ■ l' U 7 

The three eigenvectors ef^ in (|T0[) form an orthonormal set. Due to the tracelessness of S^ 2 \ the distortions associated 
with the three directions add up to zero, = —(S^ + S^). In terms of the Euler angles <E>( 2 \ @( 2 ) and \E'( 2 ), for 

(2) 

rotation of unit vectors along coordinate axes to {e^ }j=i.2,3> the eigenvectors may be represented as 

ei = [cos ^> cos $ — cos sin sin <I>, cos *Sf sin $ + cos cos ^ sin $, sin "J sin 0] , 

e2 = [— sin ^ cos $ — cos sin "J cos $, — sin \& sin $ + cos cos "J cos $, cos ^> sin 0] , 

ez = [sin© sin$, — sin© cos$, cos 0] , (71) 

where the (2) superscripts are suppressed. Just as the distortion values, the angles may depend on r. If the angles do 
not depend on r, then those angles coincide with those for the diagonalization of the tensor of quadratic moments, 
(r Ql r a2 ), and for the diagonalization of the quadrupole coefficients of the correlation function at any q. The 5 
independent parameters at any relative variable correspond to 5 independent expansion coefficients for rank I — 2. 

With regard to the higher, I > 2, shape distortions, one way to associate the directions of axes with them, is to 

it) 

choose one direction, defined by e\ , by demanding that the distortion along that direction, equal to 

s^m= £ $w.^ e w ... e w , (72) 

is extremal, minimal or maximal [25| . Directing a coordinate axis along e\ , makes the two 5^-coefficients, charac- 

(£) 

terized by £± = £ — 1, vanish. The next axis, defined by e 2 , may be chosen again as one extremizing the distortion, 
but now within the direction plane perpendicular to . When directing the second coordinate axis along e% and 
the third along the direction perpendicular to the first two, we find that the choice puts to zero the -coefficient 
characterized by £\ = 0, I2 = £ ~ 1 and £3 = 1. 
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B. Harmonic Decomposition vs Fits in Three-Dimensions 

For a given total momentum, the correlation function contains three-dimensional information pertinent to the 
three-dimensional structure of the source function. Effectively, there are three competing approaches for handling 
the information in the correlation function. In the first approach, most common till now in literature, the data is 
presented along chosen discrete directions of relative momentum, typically one or three. In the second approach 
the data is fitted to a parameterized form, e.g., a three-dimensional Gaussian, in the full three-dimensional space of 
relative momentum. When going beyond the simplest parameterizations in this approach, it may be difficult to assess, 
working in the three-dimensional space, which aspects of the source parameterization affect which particular aspects of 
the correlation. In the third approach, the data is expanded in harmonic functions. A set of one-dimensional functions 
then yields a complete representation of the three-dimensional data. Upon decomposition, the deduction of source 
features is reduced to a set of one-dimensional problems, with one-to-one correspondence between the coefficients of 
the correlation and the source. A specific coefficient of the source affects the specific coefficient of the correlation 
function and no other independent coefficients. Whether one were to advance the source analysis as an inversion 
problem, or as a problem of fitting, it should be certainly preferable to divide and conquer, i.e. approach the issue as 
a set one-dimensional problems. 

Spherical harmonic decompositions generally carry more information than can be expressed with parameters of 
gaussian sources employed in the literature. The integrated nine independent source coefficients, S^', for t < 2, fully 
determine the moments (r ai r a2 ) and (r a ), for all a-indices, cf. Eq. (f52|) . Those moments and the integral of S^°\ 
A = 4-7T f drr 2 6>(°'(r), determine all the ten parameters that can be set for a gaussian source such as discussed in [25[ 
or in the next section. However, the coefficients for £ < 2, may have a different dependence on r than those 
for a gaussian source. In addition, the £ > 2 coefficients may further yield different characteristic moments. When 
coordinate axes are directed along the principal axes of the gaussian quadrupole deformation, the gaussian moments, 
relative to gaussian center, factorize: 

((n - fa))* - (r 2 » f2 (r 3 - (r 3 »^) = <(n - (n))*) <(r 2 - (r 2 »^) <(r 3 - (r 3 »^) . (73) 

For distinguishable particles in the midrapidity region of symmetric ulrarelativistic collisions, the above implies that 
(''out r iong) = ( r out) (f\ ons ) • The latter factorization will be violated for instance for boost-invariant particle emission, 
such as discussed in Sec. IV Bl In the case of boost-invariant emission, with a transverse expansion, the single-particle 
sources acquire a boomerang-type shape, as particles at lower positions along the beam axis get emitted locally earlier 
than particles at higher positions. The violation of the factorization above will become apparent if one e.g. observes 
that the boomerang shape will be preserved in the relative source when one of the emitted particles is heavy (baryon) 
and the other light (n meson). Features of a boomerang shape may be quantified with the source expansion coefficients 
for £ = 3, which contribute to the moment (r ou t rf ong ) , cf. (|52)l . that may be compared to the product (r out ) (rf ong ) 
from I = 1 and £ = 2, besides £ = 0. 

With regard to the dependence of <S^ on r, already several studies have focused on determining the non-gaussian 
behavior of the angle-averaged source function, S<- i= V(r). With regard to £ > 1, significant values of S^(r) are 
principally expected only there where S^ ) (r) begins to fall, i.e. around the edges of the source, whereas, at small r, 
tx r . Concerning the low-r region, in the source-correlation relation, the convolution with kernel JCi brings in 
a factor r 2 . That low-r suppression generally reduces the impact of small-r on the correlation, enhancing errors 
in the restored source, making it likely difficult to detect any nongaussian features there at £ > 1. On the other hand, 
in the source tails at £ > 1 in the source coefficients, the non-gaussian behavior may be actually enhanced compared 
to £ = 0, both at low and high reaction energies. For instance, long-lived slowly cooling systems with emission rates 
falling off as e _t / r should give rise to relative sources with an exponential fall-off along the direction that the pair total 
momentum has in the frame of the emitting system. For low £, strength of the exponential tail should increase with I 
for distortion along the direction of the total momentum, Eq. (|72p . In fact, looking for the maximal distortion in 
source tail may help to locate the emitting source within the velocity space. In ultrarelativistic collisions, the emission 
becomes boost-invariant around the midrapidity. For a freezing-out boost-invariant thermal source, the single-particle 
distribution of emission points in the space-time, cf. Eq. ([2]), has the form 

s(p,r, i)~ 8(T-Tf)e- E± V^+^/T, (74) 

where r = \Jt 2 — z 1 and z points along the beam axis. The single-particle source above falls off exponentially with 
an increase in \z\ and so would the relative source with an increase in |n ng|- This behavior should be particularly 
revealed in the S[o,o,2fc] coefficients of the relative source. 

While it may be straightforward to compute the expansion coefficients as convolutions of the correlation with 
different harmonics, some caution should be exerted in selecting the range of relative momentum q for the source 
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investigation. The caution is required due to the fact that there exist other causes for inter-particle correlation than 
the final-state effects within a two-particle system. At the mathematical level, those causes can give rise to a q- 
dependence of the source S, assumed to be negligible in arriving at Eq. {]]). The dependence might be produced by 
a collective expansion of strength comparable or larger than the local thermal motion in the emitting region. If the 
dependence is not explicitly corrected for, the analysis for source restoration should be limited to the region where the 
source g-dependcnce is deemed weak. In case of the final-state effects due to resonances, short-distances in S are of 
importance, with the g-dependence governed there by local temperature at freeze-out, giving rise to the requirement 
of q <C \/2jTT for the source analysis. This condition is practically always satisfied for pp pairs whose correlation 
function peaks at q ~ 20 MeV/c. On the other hand, for pir + pairs in the region of A-resonance at q ~ 225 MeV/c, 
the condition is practically never satisfied in nuclear reactions. Correlations due to directed flow associated with the 
reaction plane, such as due to elliptic and sideward flows, can be partially accounted for by controlling the reaction 
plane in constructing the denominator in (fT]). Otherwise, those type of correlations tend to peak at q of the order 
of the mean p± (i.e. a few hundred MeV/c for p or ir in energetic collisions), reaching there values of the order of 
few percent. By comparison, Coulomb correlations fall off with increase of q, as l/q 2 , down to about one percent 
for pp typically at q ~ 100 MeV/c. That further underscores the need to limit the range of q in source analysis. In 
ultrarelativistic collisions, jets are an additional cause for inter-particle correlations. The characteristic scale for those 
correlations is of the order of several hundred MeV/c and the magnitude of the order of one percent towards lower p±. 
Overall, it is apparent that caution should be executed when considering small values of the 1Z, i.e. of the order of 
half a percent or less. If conclusions are to hinge on such small values, some supporting evidence should be provided 
to justify the assumption that the correlation originates from final-state interactions or symmetrization. 



V. EXAMPLES OF EMISSION 



Here, we provide examples of what may be expected in reactions regarding harmonic characterization of emission 
sources and of correlation functions. 



A. Gaussian Sources 



In the analyses of correlation functions, it is quite common to aim at a gaussian representation of the emission 
sources. When cast in the gaussian form, a source is expressed in terms of a minimal number of parameters needed 
for shape description. 

Thus, a gaussian source is generally described in terms of 9 parameters, out which 3 describe the position of the 
source center d and 6 correspond to the independent independent matrix elements of the symmetric positive shape 
matrix M: 

[det Ivl] | l 1 

^ ' I ai a 2 ) 

A tenth parameter would describe the normalization of the source, which could vary from unity if a significant fraction 
of the pairs had effectively infinite separations due to weak decays and were outside the description in terms of a 
single gaussian. For identical particles, the accessed relative source must be symmetric under inversion, representing 
particle interchange. In that case d = 0, leaving only 6 independent shape parameters for the source. 

The convenient choice for the shape parameters are three Euler angles describing an orientation of the eigenvec- 
tors {ui} of M and three gaussian radii Ri in terms of which the matrix is M ai a2 = ^ i (l/i?|)u ictl Ui a2 . The 
specific factor of (1/4) multiplying M and (1/i? 2 ) in (Tf5|) is associated with the tradition of assigning the radii to the 
single-particle sources s in ©. However, absolute positions in a reaction are not observable; for the description of 
relative positions more appropriate would be a factor of (1/2) yielding radii larger by a factor of V2. When coordinate 
axes are oriented along the gaussian shape eigenvectors, the gaussian shape takes the form: 

*W - (47r) 3/2 Rx Ry Rz eX P \ 4R 2 4R 2 4R 2 ]■ W 

The displacement d of the source for distinguishable particles, relative to r = 0, is of interest in both intermediate- 
energy and high-energy collisions [33l . l34j ] . For instance, at low energy it might be of interest whether neutron emission 
preceded proton emission (35| . producing a d in the direction of pair total momentum, or whether intermediate-mass 
fragments were emitted simultaneously with the protons [32j . At high energy, one might look for evidence of strange 
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particles leaving early j36j . A simultaneous emission combined with radial expansion should lead to heavy particle 
emission shifted out in the direction of expansion relative to light particle emission [3^ |. 

In many circumstances, the structure of anisotropies for measured correlation functions and corresponding structure 
for imaged sources simplifies. Thus, if the reaction plane is not identified, the measured correlation function and the 
corresponding source must obey reflection symmetry with respect to the plane formed by the pair total momentum 
P and the beam axis. In the gaussian representation, the displacement vector d must lie within that plane and so 
must two of the eigenvectors of the gaussian shape matrix M. In the midrapidity region of a symmetric system, the 
correlation function and source must be moreover symmetric with respect to forward-backward reflection. That makes 
the d- vector perpendicular to the beam axis and makes the in-plane gaussian-shape eigenvectors pointing along and 
perpendicular to the beam axis [3[. 

The nine gaussian parameters discussed above may be related both to the angular moments of the source and of 
the correlation function. It may be tempting, in particular, to associate the three independent £ = 1 moments with 
the displacement d and the five £ — 2 moments with the shape anisotropy of a gaussian. However, the situation is 
more involved. For illustration, one can consider an isotropic gaussian displaced in the x-direction: 

1 f - d) 2 + y 2 + z 1 \ 1 f r 2 +d 2 \ f xd \ 

1 f r 2 + d 2 \ ^ {2£+l)\\ ( rd \ M) fn . 

- ^j^-H^^/S^r^ 7 ^^^"^- (77) 

Here, I( are the modified spherical Bessel functions of the first kind. For small argument values of those functions, 
the functions are proportional to the power of the argument, Ii(x) oc x . Equation (|77|) illustrates that the source 
displacement d produces non-zero angular moments for all € > 1 even when the gaussian source is isotropic around 
its center. Incidentally, in a similar manner, for d = 0, the anisotropy of a gaussian gives rise to non-zero angular 
moments for all even £ > 2, not just £ — 2. 

For a general displaced and anisotropic gaussian, the angular coefficients can be obtained, as a function of r, through 
a direct angular integration. Otherwise, the power series for an exponential may be employed, in combination with 
the properties of spherical harmonics, to yield converging series for the coefficients. Specifically, as the argument of a 
source gaussian is quadratic in coordinates, that argument can be represented in terms of harmonics of rank £ < 2, 
leading to: 

5(r) = c(r) e r(r) , (78) 

where 

r(r) = T *e™(r) ^(fi) = £ 7 W h l > h 1 ' . (79) 

1<£<2 I 
m 1<^<2 

The exponential e r may be next expanded in the power series yielding, 

oo 

S(r)=Y / S N (r), (80) 

jv=o 

where 

S N (r)=S N ^(r)^, (81) 

and So(r) — c(r). Individual terms of the sum in ([80| may be represented in terms of the surface spherical harmonics 
of either type, leading to corresponding expansion for the source S. The rules of superposition for tesseral harmonics 
yield the following recursion relation for the tesseral expansion coefficients: 



Smm(r) = -j^- Y ^ ' m ' ^" ' ( m ~ m ')\ im ) r e' ™'( r ) s (N~i)e>> ( m - m >)(r) , (82) 
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where (-|-) represents a Clebsch-Gordan coefficient. It follows that terms in the sum (|80| decrease rapidly with TV, 
for iV 3> max^ m |r^ m (r)|. When aiming at the source in terms of cartesian harmonic coefficients, first the cartesian 
coefficients S N g which do not meet the tracelessness condition are found through recursion: 

= JV^5Z 7 ^ 7 ^~"^ r '*^ V-l)(/-m)( r )' ( 83 ) 



In the next stage, the traceless coefficients S^are found following such a procedure as in Eqs. (|65]) - (|66p . Finally, the 
coefficients for the source may be obtained by summing over the respective series: 

oo oc 

Se m (r) = ^2 S Ntm(r) , or Stfr) = ^ S Nii r ) > ( 84 ) 

N=0 N=0 

depending whether the tesseral or cartesian harmonics are followed. Coefficients for gaussian sources in this Section 
have been obtained following the recursion outlined above. Kernels for use in Eqs. (|4|) and (| 10[) have been obtained 
using [29] known phase-shifts. 

In order to illustrate the possibility of measuring source shapes, we consider a relative source represented by an 
axially symmetric gaussian oriented along the beam axis with parameters R x — R y = R±_ — 4 fm and R z = 8 fm. For 
distinguishable particles, we additionally displace this gaussian away from the origin at r = by d = (0, 0, 4 fm). Such 
a source may result in consequence of emission of the individual species within two single-particle sources differently 
elongated along the z axis, with the source for the a species displaced by d along the z-axis relative to the source for the 
b species. For the chosen relative source, the anisotropies of associated correlation functions can be quite significant. 
Examples of the associated correlations, for different particle pairs, pK + , p^ + , pp, pn, nn and pA, are presented in 
Fig. El The correlations are shown there as a function of cos9 q at fixed q. The observed angular dependence of the 
functions is exclusively due to the source anisotropy. While the correlations are stronger at low than at high q, the 
correlations and their anisotropies are generally of the same order. In the dependence of correlations on cos8 q , for 
distinguishable particles, one recognize both components odd and even in cos9 q . The odd components are associated 
with the displacement d in the source while even are associated with both the displacement d and R± =^ R z . In the 
correlations of identical particles, only the effects of R± ^ R z are present. 

With the choice of axial symmetry for the specific source, all cartesian components with odd number or x or y 
indices vanish. For identical particles, this further implies vanishing of cartesian components with odd number of 
z indices, since there are no indices left to pair them with to yield a net even index number. For distinguishable 
particles, at £ < 3, the only nonvanishing moments are: S^°\ <S* , Six, Syy , <S*« , S x %, S yy \ and Sj^ x . Not all of 
those moments are independent, as 5^ = — (Syy + si 2 J^j and S x 3 Jz = — (s yv \ + si 3 J z ^j , in consequence of the tensor 

tracelessness. In addition, the axial symmetry implies equality between components: S X 2 J = S y 2 y and S X xz = S yy \. 
The left panels of Fig. [5] show the £ < 3 cartesian coefficients of correlations for proton-charged meson pairs. Those 

(2) (2) 

coefficients which are not shown are easily generated from those which are, as e.g. S xx — —S zz /2, etc. 

The correlations of protons with charged mesons in Fig. [5] are plotted for q < 80MeV/c. Within that momentum 
range, the correlation should be dominated by Coulomb effects. For comparison purposes, also shown in the figure 
are the p-K + correlations expected in the classical limit of the Coulomb interactions. As expected following the 
discussions in Sec. [U the quantum effects seriously affect the correlations at qR < I + i, dampening out in that 
region the classical structures and hampering thus practical extraction of source shape information. With a smaller 
reduced mass in the jm + than in the pK + system, the classical Coulomb hole has a smaller g-size in the pn + system. In 
consequence, the quantal effects have a stronger impact on the features of the pir + correlations, making the extraction 
of source shape harder from the pm + than from the pK + correlations. 

Although difficult, it is certainly tenable to measure correlations of the order of one percent or less. In order to 
exhibit better the correlations at large q, Fig. [7] displays the same proton-charged meson correlations as those shown 
in Fig. [51 but now weighted by q 2 . It is apparent that the correlations decrease as I increases. An important feature 
of the exhibited correlations is that they are smooth in q for all I. Due to this smooth behavior, it is possible to 
analyze those correlations utilizing wide momentum bins. The phase space increasing as q 2 should further help to 
garner statistics needed for analyzing the correlations at the level of one percent or less. The practical limitation for 
the analysis at large q may turn out to be, though, the need to separate the correlations due to final-state effects 
within the pair from the possible competing correlations induced by collective flow and, at RHIC, by jets. 

Correlations of elementary baryons are good candidates for source analyses, due to the large scattering lengths within 
the pairs which lead to strong correlation structures on the relative- momentum scale of q ~ 20 MeV/c. The right-side 
panels of Fig. [5] display the correlations for pn, nn, pp, pA pairs, associated with the gaussian sources. When present, 
the £ > 1 baryon-baryon correlations tend to be significantly stronger at most momenta of experimental interest than 
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FIG. 5: (Color online) Correlations associated with a relative source in the form of an axially symmetric gaussian, characterized 
by R x = R y — 4fm and R z = 8fm, for different indicated particle pairs. For distinguishable particle pairs, the center of this 
gaussian is shifted by d = (0,0, 4 fm) relative to the origin at r = 0. The correlations are shown as a function of cosine of 
the angle that the relative momentum q makes with the z-axis, for two values of momentum magnitude, q — 25 MeV and 
q = 55 MeV. For the specific employed source, angular anisotropies in the correlation functions are of the same order as the 
functions. 



the corresponding Coulomb-dominated meson-baryon correlations. The odd correlation moments vanish for the pp 
and nn pairs as the pair sources and correlations are symmetric under particle interchange. The odd-i' correlations 
are especially important for the pn pairs emitted from intermediate-energy heavy ion collisions as these correlations 
tie to the issue of whether neutrons are emitted earlier than protons and, further, to the symmetry energy. Within 
high-energy collisions, odd-£ pA correlations are of interest as they can provide information on whether the hyperons 
leave earlier than protons. 



B. Boost-Invariant Blast- Wave Sources 



Results of early pion-pion correlation measurements at RHIC have been satisfactorily described with fits within 
the blast wave models [U [38| . Those fits have suggested a nearly instantaneous emission directly contradicting the 
formation of a long-lived mixed phase which would have resulted in a dilatory emission and in relative sources with 
sizes extended in the outward direction of the pair total momentum. Given the importance of those measurements 
and findings, it can be important to verify the conclusions on the nature of particle emission at RHIC by extracting 
sources from the correlations of other pairs such as involving kaons or protons. 

The minimal blast-wave model used here for illustration employs four parameters describing the single-particle 



24 



0.05 
-0.00 



-0.0 
-0.1 



cr 



-0.0 



0.1 - 




o - 



p-n + □ pK + - 

■ pK + (classical) _ 




R(0) _| 
i i i i I i i i i I i i i i I i i i i 



0.02 
0.00 
-0.02 



0.05 
-0.00 
-0.05 
1 





0.2 







Z e. A 


0.1 


— A a 

~ A n " 




- o ° □ a a ei 

A DO 


-0.0 


□ a o 



: „ R (3) 



ZZZ o ° ° ^ 

O 



□ o 











+H-H- 


+H-T- 


-H-H- 





2 2 2 2 ? 9JA 



— o „ z — 

o 



□ □ 



: 



□ □ □ □-! 



j i i r | 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



o 

IB □ , 



a pp □ pA — 

« nn o pn 



• £8gg gBB8 HBB-B 



R° 



a| I I I I I I I I I I I I I I I I I I 



20 40 60 80 20 40 60 80 

q (MeV/c) 

FIG. 6: (Color online) Angular coefficients of the correlations associated with an axially symmetric gaussian source characterized 
by R x = Ry = 4 fm and R z = 8 fm, for proton-charged meson pairs (symbols in the left panels) and baryon-baryon pairs (symbols 
in the right panels), as a function of the relative momentum q. For distinguishable particle pairs, the center of the gaussian 
source is shifted by d = (0,0,4fm) relative to the origin at r = 0. The pK + and piv + correlations in the left panels tend to be 
dominated by the Coulomb interaction. To illustrate the degree to which quantal effects wash out the information on emission 
source, the left panels show also (solid lines) the pK + correlations in the limit of classical Coulomb interactions. 



sources that are boost-invariant along the beam axis: 

s(p,r,t) oc 6(T-Tf)6(r max -r±)m± cosh (y-r))exp f \ . (85) 

Here, y and r\ are, respectively, energy-momentum and space-time longitudinal rapidities. The emission is assumed 
to occur at one proper time r/ = lOfm/c. The decoupling is assumed to be spread out over the cylindrical volume 
of radius r max = 13 fm and take place in the presence of a collective motion characterized by transverse rapidity p 
(u± = sin p) proportional to transverse position, p oc r±, producing the maximal transverse velocity in the longitudinal 
center of mass of w™ ax = tanh p max = 0.75 c, at r± = r max . The breakup temperature is assumed to be T — 105 MeV/c. 
The chosen parameters describe [37j particle spectra at p± < 1 GeV/c, as well as basic aspects of pion-pion correlations, 
in Au + Au collisions at ^/snn — 130 GeV. Here, our goal is not to fit data, but to use the blast-wave model to see 
what results might be expected for correlations of other particle pairs analyzed in terms of TZ^ (q) . 

We focus on correlations of nonidentical particles, because of the possibility for associated finite odd angular 
moments. The I < 3 relative-source moments for the pK and pn pairs moving at transverse velocity of v±_ = 0.6 c 
are shown in Fig. [5J To obtain the displayed results, emission positions r a j, for each of the particles a and b had 
been generated within the particle rest frame following the blast-wave description. After generating several thousands 
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FIG. 7: (Color online) The same proton-charged meson correlation coefficients as in Fig. [6] but now multiplied by the factor of 
q 2 to better exhibit the features of the correlations at large q. 



points, the number of relative positions r a — for the relative source becomes of the order of a few millions, sufficient 
for finding quite smooth source moments, except at the lowest \r a — r&|. We associate the z-axis direction with the 
beam axis and the x-axis direction with the outward direction along the pair momentum. The angular moments for 
the blast- wave model are compared in Fig. [8J to the moments for a Gaussian source with the same cartesian moments, 
(x), (x 2 ), (y 2 ) and (z 2 ), as the blast- wave relative source. 

As elsewhere, values for the underlying three-dimensional source structure can be restored along any direction from 
the cartesian coefficients, in particular from those displayed in Fig. [5J exploiting tracelessness. Thus e.g. along the 
x-axis the source values are 

S = ± + — \S^J — S^J^J — 2^zz "F &xyy "F S^hzz + • • • > (86) 

where the upper and lower signs refer to the positive and negative directions along the axis. Within the xy-plane, at 
45° to the x-axis, the source values are 

S = ± ± SP - \ Sf) T ^ Sit + -., (87) 

with the upper signs again referring to the positive direction along the x-axis. 

Within the blast-wave model, the heavier particles end up being emitted farther out, on the average, along the 
outward direction, than the lighter particles. This is due to the following. In an infinite system with Hubble expansion, 
a particle with a given velocity would get, on the average, emitted from a spatial patch centered around the point 
with the collective velocity matching the particle velocity. The size of such an emission patch would depend on the 
particle mass, with the patch being tighter for heavier particles with lower thermal velocities and being more extended 
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FIG. 8: (Color online) Angular cartesian coefficients, at £ < 3, multiplied by 47rr 2 , of relative sources for pK (left) and pir 
(right) pairs moving out at transverse velocity of 0.6 c relative to the beam axis. Symbols represent the results obtained in 
the blast-wave model. Lines represent the results from a gaussian source with the same pair cm. cartesian moments, (n) and 
(ri Vj), as the blast- wave source. The non-gaussian features of a blast- wave source become more apparent as the mass difference 
within a pair increases. 



for lighter particles. Relative to the infinite system, however, in a system with a finite transverse extension, the outer 
portions of the emission patches would get cut off and more so for the light than for heavy particles, forcing the 
emission region of the light particles inwards relative to the heavy. In the past, relative displacements of emission 
between particles of different mass, e.g. Ktt, have been accessed by comparing the correlations at opposite signs of 
outward momentum, q out > and g out < 0, at low values of the remaining q-components [HI, |4(j. Such a strategy 
discards potentially useful information on the relative displacement in the remainder of the relative momentum space. 
The displacement of species within a pair is visible in Fig. [8] in the finite values of . With the mass difference 
being larger within the pir than pK pairs, the coefficients peak farther out and reach significantly higher values 
for the pir than pK pairs. 

Some other features of blast-wave sources become apparent in the context of comparisons to gaussian sources. 
When making such a comparison, a few points need to be kept in mind. In an infinite system with a Hubble-like 
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expansion along the longitudinal direction, the associated relative source should fall off exponentially at large r in 
the longitudinal direction, rather than falling in a gaussian fashion, see e.g. Eq. (j?4")) . On the other hand, when the 
emission within a model gets confined to a finite transverse domain, the source function should fall off more sharply 
in the transverse directions, especially in the sideward direction which is less affected by decays. For our comparison 
we choose the parameters of the gaussian source to match the norm and (r^) and (r^ rj) of the blast- wave source. 
Given Eq. (|53|) , this means that the £ < 2 integrals of r e+2 S^(r) for the gaussian are made to match those for the 
blast- wave. 

When comparing the blast-wave and gaussian source-coefficients in Fig. [5J we can see that an agreement between 
them generally worsens as £ increases and it is worse for pir than pK pairs. These tendencies are to an extent 
associated with the source tails. Since the emission region region is wider in the case of pions than kaons, the larger- 
scale geometry of the emitting region is explored more stringently in the case of pions. In addition, the pions become 
relativistic at lower pt than do kaons. The finite transverse extension of the emission zone most strongly suppresses 
the tail of the relative source in the negative x-direction. In terms of angular coefficients, this suppression generally 
results in a large enhancement of the high-£ blast-wave coefficients at high-r, compared to the gaussian coefficients. 
This enhancement is particularly visible for the pn pairs. The one case in Fig. [5] where the pn S zz is smaller in 
magnitude for the blast-wave than gaussian source is associated with the change in sign for that coefficient. The 
change in sign is itself associated with the trimming of the relative source by extension of the emitting region, present 
in the x and y directions, but absent in the z-direction. Further sign changes are seen for the blast-wave £ — 3 
coefficients, while there are no such changes for the corresponding gaussian coefficients. These sign changes reflect the 
more complex nature of the blast-wave sources of a crescent-like shape in the xy and xz directions. The latter shape 
results from the relativistically noninstantaneous nature of the boost-invariant blast-wave emission, mentioned already 
in the preceding section. Coded in the Sx yy and S xzz coefficients, cf. (|53p . are the violations of the factorizations which 
hold for the gaussian source: (xy 2 ) — (x) (y 2 ) and (x z 2 ) = (x) (z 2 ). These violations result in significant differences 
between the £ = 3 coefficients for the gaussian and blast wave. 

Regarding the examples of source coefficients in Fig.[5J it is important to note that, while the £ = coefficients have 
most of their strength within the (5 — 15) fm region of separation, the £ > 1 coefficients have most of their strength 
within the (7 — 25) fm region. The fact that the source anisotropies are most pronounced at large r, combined with 
the difference in fall-off for strong and Coulomb kernels, of 1/r 2 vs 1/r, allows the Coulomb induced correlations to 
compete effectively with the correlations induced by strong interactions in providing access to the source anisotropies. 

Figure [5] displays the proton-charged meson correlations calculated by convoluting the source functions illustrated 
in Fig. [8J with the respective kernels. The utilized pK + and pn + kernels have been obtained accounting both for the 
Coulomb as well as the strong interactions [2^|. Relative magnitudes of discrepancies between the blast- wave and 
gaussian are similar for the corresponding correlation and source coefficients. The level of difficulty in accessing the 
shape information may be assessed in Fig. [9] by examining the absolute magnitudes of coefficients for the different £- 
values. For I = 1, the magnitude is of the order of few percent, which should be observable with statistics characteristic 
for the past correlation measurements. For £ > 2, the magnitude is of the order of one percent and less. Particularly 
well suited for investigating the high-i? correlations should be the high-statistics data from the 2004 run at RHIC. 
Eventually, the reach of the analysis may be limited by systematics rather than by statistics. 

VI. SUMMARY 

Correlation measurements are a crucial tool for understanding the dynamics of central nuclear reactions. One of 
the forefront tasks before the correlation analysis is to decipher the mire of RHIC data. Without carrying out a 
multidimensional correlation analysis, there is little chance for a detailed understanding of the impact of dynamics on 
correlations in those reactions and for understanding of the role played by the equation of state. To date, correlation 
analyses at RHIC have focused on identical pions. In identical-pion correlations, strong interactions play little role. 
In the analyses to date, Coulomb interactions have been normally corrected for in a shape-independent manner, 
leaving off correlations induced by identical-particle interference. The latter are related to the source through Fourier 
transformation allowing for easy modelling with specific sources such as gaussian. The decomposition of correlations 
in terms of cartesian coefficients offers an alternative way of carrying out a multidimensional correlation analysis for 
identical pion pairs and it provides possibility for carrying out multidimensional analysis for other classes of final-state 
effects. 

One of the benefits of decomposing information with angular harmonics is that the correlations induced by Coulomb 
and strong interactions carry, on their own, some shape information that can be exploited in the correlation analysis. 
The representation of the correlations in terms of the cartesian or tesseral coefficients carries the full three-dimensional 
information in terms of a set of one-dimensional functions. Data represented in this fashion can be then be compared 
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FIG. 9: (Color online) Cartesian coefficients for the correlations of pK + (left) and piv + (right) pairs moving out at transverse 
velocity of 0.6 c relative to the beam axis. Symbols represent the results obtained in the blast-wave model, while lines represent 
the results from a gaussian source with the same pair cm. cartesian moments, (n) and (rirj), as the blast-wave source. The 
corresponding cartesian source coefficients have been represented in Fig. [8] As in the case of Fig. [8] differences between the 
blast-wave and gaussian results are more pronounced for the pjr + than for pK + pairs. 



to models visually, one harmonic at a time, providing the means to better identify which shape characteristics of the 
data are being described by the model. 

The ability to discern information about source shapes from correlation measurements depends on the structure 
of the kernels, }Ct(q,r), for different spherical ranks I. In Section [HI it was demonstrated that the identical-particle 
interference, Coulomb and strong interactions all produce effective kernels. The resolving power of Coulomb interac- 
tions was shown to rise with mass and charge within the interacting pair. Strong interactions are most effective in 
the resonance region and, otherwise, for large cross sections. Even purely s wave interactions can produce kernels of 
good resolving power at higher I. 

While expansions employing tesseral harmonics nominally provide as powerful means for storing information in 
correlation functions as cartesian harmonics, they do not provide as transparent or as intuitive means for understanding 
the shapes as expansions done using the cartesian harmonic basis. Cartesian harmonics have a variety of useful 
properties which can be exploited in manipulating the angular information as was shown in Sec. IIIII The connection 
between specific aspects of the source geometry and specific angular moments has been emphasized in Sec. IIVI It has 
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been pointed out that the average relative displacement of emission points in the cm. frame of nonidentical particles 
can be assessed through the £ — 1 angular moments. The quadrupolc distortion, most commonly quantified in terms 
of the radii i? ou t, ^sidc and Ri ons , can be studied with the £ = 2 moments. The £ — 3 moments test the more complex 
aspects of sources such as boomerang features. 

Detailed examples of sources and correlations have been presented in Sec.fVj Using gaussian sources, it has been, in 
particular, demonstrated that the pK + and pn + correlations can realistically be expected to provide shape information 
on emission sources. The strengths of correlations characterizing shape anisotropies can reach the magnitude of few 
percent. While the upcoming RHIC data sets will provide sufficient statistics to assess correlations between different 
particles at the level of one percent or better, the true challenge may turn out to the subtraction of competing 
correlations, such as those produced by jets. The proton-charged meson correlations are dominated by Coulomb 
interactions. The strong interactions in baryon-baryon systems can also, though, provide good means for discerning 
shape anisotropies at low relative momenta, yielding correlations of the same order or larger than the proton-meson 
correlations. One common feature of all pairs and sources examined here is that the variation of the correlation with 
angle is comparable in magnitude to the angle-integrated correlation. Consequently, if one has sufficient statistics to 
measure a one-dimensional correlation, only a modest enhancement in statistics would suffice to find, or constrain 
in a meaningful manner, the higher angular moments. In the context of the blast-wave sources, opportunities for 
detecting more complex shapes have been discussed, such as boomerang shapes. 

In conclusion, a treasure trove of potential information lies largely ignored in the huge volumes of data from high- 
energy collisions. This information can address what is the perhaps the greatest surprise of the first several years of 
RHIC experiments, the source sizes and shapes inferred from identical-pion correlations. The main difficulty with the 
proposed new class of correlation analyses is the fact that they rely on observables that are quantitatively small, of 
the order of a percent. This will require a careful analysis of competing correlations from jets and collective flow, but 
we believe that the objectives will indeed prove feasible. 



APPENDIX A: PROJECTION OPERATOR V 



In arriving at the cartesian harmonics, an important role has been played by the cartesian projection operator V 
which projects out the symmetric traceless portion of any cartesian tensor. As a projection operator, the operator V 
must have the following properties: 

v (t:t) = v {t:i) (A1) 

v (l:l) \ ' v (i:l) ^ v (t:l) , ^ 

The projection operator acting on any tensor T must yield a traceless result VT . Since we can choose, in particular, 
a tensor that is finite for any single set of indices and zero otherwise, we find that the operator V itself must be 
traceless: 

y v {i-.t) i= y V m , , =0 . (A3) 

a a 

We shall demonstrate existence of the operator V by construction. Notably, when acting onto a symmetrized traceless 
tensor, the operator must reproduce the tensor. Thus, we can represent V as 

v m) = s m) Av (t.t) (M) 

or, in the shortened notation, 

where S is the symmetrized identity operator 

S < a 1 la e :a' 1 ...a' e = ~J\ ^ ^Mi) °>'l " " " ^M«) a 't 1 (^6) 

and AT is an operator which yields a zero when acting upon a traceless tensor. In (|A6[) . the sum is over permutations 
of £ indices. Importantly, a traceless symmetric tensor of rank £' is traceless and symmetric within any subset £ < £' 
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FIG. 10: Diagrams illustrating the construction of the projection operator "p( M >, using the lower-rank operator 11 r \ 
represented by the rectangular block, and 5-symbols, represented by lines, (a) Diagrammatic representation for the leading 
term I 'p( t ' e \ Eq. (|A8[l . (b) Diagrammatic representation for the terms resulting from evaluating the trace of Ir p( e:e \ The terms 
represented by the left diagram yield zero due to the tracelessness of "p( f_1:£_1 '. (c) Diagrammatic representation for the 
n-pH-t) anc [ lil-p(l-.t) terms within the construct, Eqs. (|A9|I and (IA10|) . respectively, (d) Representation for the nonvanishing 
terms from evaluating the trace of III V ■ (e) Diagrammatic representation for the final IV/ p( e '- e > term within the construct, 
Eq. J5T1). 



of its indices. When acting on any of £ indices of a rank-£', £' > £, traceless tensor, the constructed operator p^'-^ 
will reproduce that tensor and AP^'^ alone will annihilate it. 

In our construction of V, we will find a graphical method similar to that employed for interactions [4l| useful. 
Within the method, the Kronecker 5 aa ' is represented as a line joining the indices a and a'. A convolution of delta 
symbols a <^*i a 2 ■ ■ ■ $a n a' = $aa> is represented by joining of the line segments which produces a single 

line from a to a'. A line which closes on itself represents a scalar factor which is the trace of the (5-symbol, equal to 
the dimension of space, J2 a ^ aa = 3- Examples will be provided in Fig. [TU] 

For I — and £ = 1, the operator V is simply the identity operator, p( 0:0 ) = 1 and = 5 aa >. We shall show 

that, when the operator for the rank (£—1) tensors is known, it is possible to construct the operator with the required 
properties for the next rank £ tensors. Thus, the operator for any rank can be constructed by recursion. Let us assume 
that we know the operator for rank {£ — 1), pC^- 1 ^- 1 ) . We shall construct the next rank operator as a sum of terms, 

pm = iff* + npM + ..., (A7) 

exploiting the operator p^ 1 '-^ 1 ) . The terms of different structure will be added to achieve in (|A7|) the 
properties (|A1|) - (| A3|l for p™'-*-). In particular, only I P in our construction will contain the symmetrization opera- 
tor S on the r.h.s. of Eq. (|A5|) . 

Since, within each rank, the operator must, on one hand, contain a symmetrized identity operator and, on the other, 
be traceless, it is natural to start with I P as a tensor product of the lower-rank operator and identity, illustrated in 

Fig-nnta), 

<W P ( '- 1; '- 13 , , , (A8) 
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symmetrized separately in the index sets a and a'. Evaluation of the trace for l V , produces terms which correspond to 
the trace of the lower-rank operator •p(^ _1:£_1 ) ; represented by the left diagram in Fig. fTUT b). and to the convolutions 
of the 5-symbol with represented by right diagram in Fig. [TOT b). Terms of the first type yield zero due 

to the tracelessness of 'p( e ~ 1 ' J ~ 1 \ but those of the second type generally yield a finite contribution. Aiming at the 
tracelessness of one can next add to the construct a counter-term 11 V of the form 

-c n 8 ai ia ,P ( ' _1; ' _1) > / / , (A9) 

to cancel the last contribution, symmetrized separately in the indices a and a', represented by the left diagram in 
Fig. HOT c). The constant cu should be adjusted to achieve the tracelessness of V^ e ' e \ Taking the trace of TI P produces 
the same type of terms as for 11 V, that were represented in Fig. llOf b). Thus, an adjustment of c/j can indeed produce 
a traceless sum T V + II V. However, the operator V^ 1 '-^ should be symmetric under the transposition of the indices 
a and a', Eq. (jAip . In consequence, when adding the term 11 V to the construction, one also needs to add a term 
IH-p = (U-p) 1 1 of the form 
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-l;t—l) 
ai...at- 1 :a' 1 ...a' e _ 2 at 



(A10) 



symmetrized separately in the index sets a and a' and represented by the right diagram in Fig. llOf c). The nonvanishing 
terms from evaluating the trace of 111 V arc illustrated with the diagram in Fig. [TOT d). That last contribution to the 
trace may be compensated for by adding the final term IV V to the construct, symmetric under transposition, of the 
form 

CIV Sae-! at ^a' e i a' e ^ ^'ai...o/_2 0-a' 1 ...a' l _ 2 ' (AH) 



symmetrized separately in the indices a and a', and illustrated with the diagram in Fig. HOf d). Taking the trace of 
1V V produces, up to a factor, the same result as taking the trace of III V . 

The constants en and cry can be now set to make the trace of "P^W) vanish. The net result for the operator, with 
the symmetrizations made explicit, is 



■pirn 



a 1 ...a t :a , 1 ...a , ll g 

m,n— 1 
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p2 / j a m a n ai ...a m -! a m+1 ...a e :a' 1 ...a' n _ 1 a' n + 1 ...a' e 



£ 2 {21 - 
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P- {21 - 
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/ j e*m oir cKi-..a m _i CK m _)_i...a, — x cK r _|_i...a^ a f n :a , 1 ...a , n _ 1 at' n , 1 ...atg 
<m<r<l 
l<n<£ 



Km<e 



+ £ 2 (2e-l) 2 ^ S <x™a r $a> n a> g 
^ ' l<m<r<e 

l<n<s<£ 

y (t-iu-i) (A12) 

a 

By now, we have constructed an operator "P™;£) wn j cn j s symmetric separately in the covariant and contravariant 
indices, traceless and symmetric under the transposition. The remaining question is whether repeated applications of 
the operator yield the same result, i.e. whether Eq. (|Al|) is satisfied. 

To answer the last question, let us apply J>( e '< e \ out of the four components in Eqs. (|A7|) and (|A12[) . to the constructed 
r p( e '' e '> with the properties we have established. The I V component reproduces P^', as both the Kronecker symbol 
and "P™ -1 ^ -1 ) reproduce a traceless symmetric tensor. On the other hand, all the other components of the construct 
annihilate V^' ,l \ because they involve calculation of the trace of that operator, either directly or after application of 

Use of Eq. (|A12|) for £ = 2 produces 

^ii a 2 -a[ a' 2 ~ o (^ Ql a 'l a 2 Q 2 ^ Q 2 a 'i) ~ o a 2 a' 2 ■ (A-13) 
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In general, for 'p^'^ at I > 2, the repeated recursion with Eq. (|A12[) produces, besides the leading symmetrized 
identity term, the correction terms which involve replacing k pairs of Kronecker symbols in the identity, linking 
the covariant and contravariant indices, with the pairs of Kronecker symbols linking separately the covariant and 
contravariant indices, k = 1,2 .. .1/2. Application of both sides of Eq. (|A12[) for V to n e produces the recursion 
relation ((48]) for A. 
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